Historic,  Archive  Document 

Do  not  assume  content  reflects  current 
scientific  knowledge,  policies,  or  practices. 


aS494 

.5 

. E8A8 

* r / i u 


TECHNICAL  MEMORANDUM 
NO.  9 


AN  EVALUATION  OF  RELATIONSHIPS  BETWEEN  VEGETATIVE  INDIES, 
SOIL  MOISTURE  AND  WHEAT  YIELDS 


UNITED  STATES  DEPARTMENT  OF  AGRICULTURE 
FOREIGN  AGRICULTURAL  SERVICE 

CROP  CONDITION  ASSESSF€NT  DIVISION 
HOUSTON,  TEXAS 


0irp  I 8 1979 


NASA-JSC 


TM-9 


UNITED  STATES  DEPARTMENT  OF  AGRICULTURE 
FOREIGN  AGRICULTURE  SERVICE 


AN  EVALUATION  OF  RELATIONSHIPS 

BETWEEN  VEGETATIVE  INDICES,  SOIL  MDISTURE 

FIRST  ISSUE 

AND  WHEAT  YIELDS 

Approved  By: 


Jirrmy  D.  Murphy 
Director,  Crop  Co 


1.  REASON  FOR  ISSUANCE 

Document  the  results  of  this  study  for  use  by  the  Crop  Condition 
Assessment  Division  (GCAD)  of  the  Foreign  Agricultural  Service 
(FAS).  The  results  presented  in  this  dociment  serve  as  a 
foundation  to  future  research  aimed  at  identifying  predictive 
relationships  between  Landsat  rrult  i spectral  scanner  data  and 
ground  observed  measurements  such  as  plant  density,  plant  height, 
and  yield. 


2.  COVERAGE 


The  paper  presents  the  relationships  among  several  vegetative 
indices,  relationships  between  the  vegetative  indices,  soil 
moisture  and  a number  of  ground  observed  data  including  plant 
density,  plant  height  and  yield  of  spring  and  winter  wheat. 
These  relationships  were  derived  with  respect  to  the  growth 
development  curve  of  vheat. 
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PART  1.0  INTRODUCTION 

1.1  SIM4ARY  AND  CONCLUSIONS 

1.1.1  Surrmary . The  mission  of  the  Crop  Condition  Assessment  Division 
(OCAD)  of  the  Foreign  Agricultural  Service  (FAS)  is  to  verify 
and  assess  the  impact  of  adverse  conditions  on  crops  in  impor- 
tant foreign  producing  areas  and  report  the  results  in  a timely 
manner  to  FAS,  Washington,  D.C.  CCAD  crop  condition  assessment 
reports  are  used  by  FAS  ccrrmodi  ty  analysts,  along  with  other 
information  sources  to  shape  U.S.  foreign  trade.  The  timely 
reporting  of  crop  condition  information  can  be  translated  into 
substantial  benefits  to  U.S.  foreign  trade,  the  American  farmer 
and  consumer. 

The  CGAD  uses  many  different  types  of  data  sources  to  support 
its  analysis  including  Landsat,  meteorological,  soils,  and 
other  ancillary  data  sources. 

The  CCAD  uses  Landsat  MSS  digital  data  in  determining  the  areal 
extent  of  an  abnormal  event  affecting  agricultural  production 
and  in  assessing  crop  condition.  Vegetative  Indices  (Vi's)  are 
computed  frcm  the  MSS  digital  data  and  are  used  by  OCAD  analysts 
to  measure  the  relative  "greenness"  of  agricultural  areas.  Re- 
searchers have  developed  a number  of  Vi's,  all  of  vhich  are 
sensitive  to  the  density  of  the  vegetative  ground  cover.  The 
basis  for  development  of  these  indices  was  to  reduce  the  four 
Landsat  spectral  bands  to  their  greenness  component . The  greater 
the  density,  or  canopy  of  vegetative  areas,  the  greater  the  VI. 
CCAD  analysts  make  qualitative  measurements  of  crop  condition 
by  comparing  the  Vi's  computed  from  different  areas  of  a region 
or  country  or  frcm  the  same  area  for  different  years.  The 
purpose  of  this  study  is  to  provide  the  OCAD  a basis  for  devel- 
oping a capability  of  using  the  Vi's  along  with  other  supportive 
data  sources  including  soil  moisture  and  crop  calendar  to  accu- 
rately measure  the  condition  or  yield  of  a crop. 

The  objectives  of  this  study  were: 

1.  To  identify  the  VI  (s)  which  provide(s)  the  most  significant 
relationship  to  plant  height,  plant  density  and  yield. 

2.  To  identify  the  time(s)  in  the  growing  season  during  which 
the  relationship  between  the  VI  (s)  and  plant  height,  plant 
density,  soil  moisture  and  yield  is  (are)  the  strongest. 

3.  To  identify  the  functional  relationship  of  the  Vi's,  soil 
moisture,  and  crop  calendar  to  plant  height,  plant  density 
and  yield. 
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The  VI 1 s*  evaluated  in  the  study  included  the: 


1) 

Ashburn  Vegetat i ve  Index 

(AVI) 

2) 

Di f ference  Vegetat ive  Index 

(DVI) 

3) 

Green  Vegetat ive  Index** 

(GVI) 

4) 

Kauth  Vegetat ive  Index** 

(KVI) 

5) 

Leaf  Area  Index 

(LAI) 

6) 

Perpendicular  Vegetative  Index  (Band  6) 

(PVI6) 

7) 

Perpendicular  Vegetative  Index  (Band  7) 

(PVI7) 

8) 

Transformed  Vegetative  Index  (Band  6) 

(TVI6) 

9) 

Trans  formed  Vegetat ive  Index  (Band  7) 

(TVI7) 

The  data  set  included  ground  observed  and  Landsat  digital  data 
collected  over  22  blind  sites  located  in  North  Dakota  and  Montana 
during  the  1978  spring  and  winter  wbeat  crop  seasons.  Ground 
observed  data  were  collected  for  10-15  \dieat  fields  per  segment 
including  wheat  type,  plant  height,  growth  stage,  yield  and  field 
ccrrments.  Landsat  data  acquired  between  April  and  Septerber  of 
1978  yielded  103  useable  acquisitions.  Image  corrections  were 
exercised  on  the  Landsat  data  prior  to  computing  the  Vi's.  These 
corrections  included:  Landsats  2 and  3 calibration,  sun  angle 
correction  to  a 39°  solar  zenith  angle,  and  haze  correction  using 
the  XSTAR  algorithm  developed  by  the  Environmental  Research 
Institute  of  Michigan  (ERIM)  . 

Soil  moisture  estimates  were  added  to  the  data  set  to  determine 
its  relationship  to  the  ground  observed  data  and  the  Vi's.  It  was 
hypothesized  that  the  soil  moisture  estimates  would  also  serve  to 
strengthen  the  relationship  between  the  ground  observed  data  and 
the  Vi's.  Meteorological  data  collected  at  several  weather 
stations  identified  near  each  of  the  22  blind  sites  were  used  to 
estimate  soil  moisture.  The  Two  Layer  Soil  Moisture  Model,  used 
by  the  National  Oceanographic  and  Atmospheric  Administration  (NCAA), 
was  used  to  compute  estimates  of  surface  and  subsurface  soil  mois- 
ture. The  model  uses  long-term  temperatures,  current  temperatures 
and  precipitation  and  potential  water-holding  capacity  of  the  local 
soils  to  estimate  soil  moisture  levels. 

Correlation  and  regression  analyses  were  performed  on  the  data  set 
to  determine  the  relationships  among  the  ground  observed  and  soil 
moisture  variables  and  the  Vi's.  The  data  set  was  subdivided  into 
seven  growth  stage  intervals  based  on  the  growth  stage  information 
obtained  for  each  field  within  each  acquisition. 


* The  VI  transformat  ions  are  shown  in  Appendix  A. 

**  The  GVI  is  equivalent  to  the  greenness  component  of  the  Kauth  trans- 
form, and  the  KVI  is  equivalent  to  the  Green  Number  used  during  the 
Large  Area  Crop  Inventory  Experiment. 
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The  growth  stage  intervals  are  shown  below. 


INTERVAL 

GRGVTH  STAGES 

LABEL 

1 

0.0 

Plant ing 

2 

1.0  - 5.0 

Til lering 

3 

6.0  - 10.0 

Stem  Extension 

4 

10.1  - 10.5 

Heading 

5 

10.51  - 10.54 

Flowering 

6 

11.1  - 11.4 

Ripening 

7 

12.0 

Harvest 

The  tillering  through  ripening  growth  stage  intervals  are  based 
on  the  Feekes  scale;  while  the  planting  and  harvest  intervals 
were  added  for  the  purposes  of  this  analysis  and  were  given  the 
0.0  and  12.0  growth  scale  definitions,  respectively. 

The  analysis  evaluated  the  relationships  among  the  ground  observed 
and  soil  moisture  variables  and  the  Vi's  at  each  of  the  seven 
growth  stage  intervals.  It  was  hypothesized  the  relationship  bet- 
ween a VI  and  crop  condition  parameters  (such  as  yield)  within  a 
growth  interval  would  more  closely  approximate  a linear  relationship. 

The  correlation  results  showed  that  the  strength  of  the  relationship 
between  yield  and  each  of  the  indices  at  each  of  the  growth  stage 
intervals  was  approximately  the  same.  The  correlations  were 
highest  and  less  variable  during  the  tillering,  stem  extension  and 
heading  growth  stage  intervals.  The  average  correlation  coeffi- 
cients between  the  Vi's  and  yield  were  -.31  at  planting,  -.26  at 
tillering,  .36  at  stem  extension,  .58  at  heading,  .17  at  flowering, 
.01  at  ripening,  and  -.42  at  harvest.  The  strongest  relationships 
between  final  wheat  yields  and  the  indices  occurred  at  heading. 

Stratification  of  the  data  by  wheat  type  (winter  or  spring  wheat) 
or  agrophysical  unit  (APU)  did  not  significantly  affect  the 
relationships  between  yield  and  the  Vi's.  Some  improvements  on 
the  relationship  occurred  at  stem  extension,  but  the  relationship 
remained  relat ively  weak.  Stratification  of  the  data  by  field 
size  did  significantly  improve  the  relationship  between  the 
indices  and  yield.  The  greatest  improvements  were  realized  at 
heading. 

Correlation  analysis  found  a strong  relationship  between  plant 
density,  plant  height  and  the  indices  at  tillering  and  stem 
extension.  Problems  with  the  ground  observed  data  prevented 
analysis  at  the  later  growth  stage  intervals. 

Subsurface  soil  moisture  showed  a stronger  relationship  to  yield 
than  was  found  with  surface  soil  moisture.  The  relationship 
peaked  at  heading  producing  a correlation  coefficient  of  .59. 

The  relationship  remained  strong  averaging  above  .45  at  the  other 
growth  stage  intervals,  except  at  stem  ext ens ion  when  it  dropped 
to  .15. 


SEP  18  1979 


3 


< 


* 


TM-9 


The  strongest  relationships  determined  by  the  regression  ainalysis 
occurred  at  heading  and  flowering.  The  coefficient  of  determi- 
nation (R^)  and  standard  deviation  between  the  observed  and 
estimated  yields  averaged  .57  and  5.75,  respectively,  at  heading 
and  .62  and  6.00,  respectively,  at  flowering.  The  relationships 
at  heading  irrproved  significantly  by  imposing  a constraint  on 
minimum  field  size.  At  heading  the  coefficient  of  determination 
and  standard  deviation  irrproved  to  .81  and  4.08,  respectively, 
after  imposing  a 30  acre  field  size  constraint.  The  relationship 
remained  stable  after  imposing  higher  minimum  acreage  limits  of  35, 
40  and  45  acres. 

1.1.2  Conclusions. 


° The  vegetative  indices  are  highly  correlated  to  each 
other  at  each  of  the  seven  growth  stage  intervals. 

° The  vegetative  indices  are  correlated  to  wheat  yield 
similarly  at  each  of  the  seven  growth  stage  intervals. 

The  degree  of  correlation  varies  by  growth  stage  interval. 

° Each  vegetative  index  is  correlated  to  plant  height  and 
density  similarly. 

° The  relationship  between  wheat  yield  and  the  vegetative 
indices,  soil  moisture  and  crop  calendar  was  strongest 
at  heading  and  flowering. 

° There  are  strong  relationships  between  plant  height  and 

density  and  the  vegetative  indices,  soil  moisture  and  crop 
calendar  at  the  tillering  and  stem  extension  growth  stage 
intervals.  The  strongest  relationships  were  found  at  the 
stem  extension  growth  stage  interval. 

° An  accurate  crop  calendar  for  wheat  is  necessary  before 
attempting  to  estimate  plant  density,  plant  height  or 
wheat  yield  using  vegetative  indices  and/or  soil  moisture. 

1.2  BACKGROUND 

The  Crop  Condition  Assessment  Division  (OCAD)  of  the  Foreign  Agricultural 
Service  (FAS)  is  responsible  for  verifying  and  assessing  the  impact  of 
adverse  conditions  on  crops  in  important  foreign  producing  areas  and 
reporting  the  results  to  FAS  ccmrodi  ty  analysts  in  Wadshington  D.C. 
Foreign  crop  areas  have  been  identified  by  FAS  that  have  a significant 
effect  on  U.S.  foreign  trade  and  are  therefore  candidate  areas  for 
timely  analysis  by  the  CCAD.  CCAD  attention  may  be  drawn  to  areas 
undergoing  abnormal  conditions  by  its  continual  screening  of  Landsat 
and/or  meteorological  data.  Additionally,  the  CCAD  may  be  alarmed  of 
events  by  information  sources  external  to  the  CCAD,  such  as  by  foreign 
agricultural  attaches  assigned  to  these  areas.  For  example,  the  CCAD 
may  be  alerted  to  a potential  drought  in  the  U.S.S.R.  vhich  could  have 
a substantial  impact  on  the  Soviet  wheat  crop  and  its  future  production. 
The  CCAD  will  provide  timely  reports  to  FAS  in  Washington  on  the 
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condition  of  the  wheat  crop,  indicating  the  drought  severity  and  areal 
extent.  FAS  carmodi ty  analysts  will  use  this  information  along  with 
other  information  sources  to  assess  the  probable  economic  impact  on  the 
world  wheat  market  relative  to  U.S.  interests. 

Methods  used  by  the  OCAD  for  assessing  crop  conditions  using  Landsat 
mul t i spectral  scanner  (M5S)  data  include  visual  screening  of  the  imagery 
and  ccnputer  analysis  of  the  digital  data.  Visual  assessments  may  be 
qualitative  statements  describing  the  condition  of  the  crop,  such  as 
the  crop  looks  better  or  worse  than  a nominal  or  reference  year.  The 
OCAD  currently  uses  NGS  digital  data  in  determining  the  areal  extent 
of  an  abnormal  event  affecting  agricultural  production  and  in  assessing 
crop  condition.  Vegetative  indices  (Vi's)  are  computed  frcm  the  KGS 
digital  data  and  are  used  by  GCAD  analysts  to  measure  the  relative 
"greenness"  of  agricultural  areas.  Researchers  have  developed  a number 
of  VI 1 s all  of  which  are  senstitive  to  the  density  of  the  vegetative 
ground  cover.  The  basis  for  development  of  the  Vi's  was  to  reduce  the 
four  Landsat  spectral  bands  to  their  greenness  component.  The  greater 
the  density  or  canopy  of  vegetative  areas  the  greater  the  VI.  OCAD 
analysts  make  qualitative  measurements  of  crop  condition  by  ccmparing 
Vi's  computed  frcm  different  areas  of  a region  or  country  or  from  the 
same  area  for  different  years.  The  purpose  of  this  study  is  to  provide 
the  OCAD  a basis  for  developing  a capability  of  using  the  Vi's  along 
with  other  supportive  data  sources  such  as  soil  moisture  and  crop 
calendar  data  to  accurately  measure  the  condition  or  yield  of  a crop. 

For  many  years,  remote  sensing  and  agricultural  specialists  have 
believed  that  predictable  relationships  exist  between  crop  vegetation 
and  the  spectral  data.  Hass  e_t  al_  (Exhibit  1,A)  used  Landsat  data  to 
qualitatively  and  quantitatively  analyze  the  amount  and  seasonal  con- 
ditions of  rangeland  vegetation.  Quantitative  measurements  of  rangeland 
condition  were  made  by  normalizing  the  difference  between  the  visible 
red  and  infrared  channels  of  Landsat.  In  essence,  Landsat  was  measuring 
the  amount  of  green  biomass  of  the  rangeland. 

Colwell  (Exhibit  1,B)  used  Landsat  data  to  identify  annual  grasses 
experiencing  moisture  stress  that  could  lead  to  potential  fire  hazards. 
In  addition,  Landsat  data  were  used  to  identify  vegetation  affected  by 
frost  damage.  The  study  objective  was  to  give  early  warning  of  condi- 
tions and  areas  that  would  require  remedial  measures  to  prevent  range 
and  forest  fires. 

The  General  Electric  Company  used  Landsat  data  to  identify  the  temporal 
spectral  profile  of  com  and  soybeans  in  Iowa  for  the  purpose  of 
investigating  crop  development,  condition  and  yield  (Exhibit  1,0). 
Promising  results  were  obtained  in  correlating  Landsat  data  to  corn 
yield.  Kansas  State  University  used  Landsat  data  to  estimate  the  leaf 
area  index  of  wheat  which  is  used  in  an  evapotranspirat ion  model. 

The  evapotranspirat ion  estimate  is  then  used  in  a model  to  predict  the 
yield  of  winter  wheat  (Exhibit  1,D).  The  results  showed  that  predicted 
wdnter  wheat  yield  agreed  favorably  with  the  average  county  yields 
estimated  by  the  Economics  and  Statistics  Cooperative  Service. 
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Thcrrpson  found  a high  correlation  between  the  Green  Number  (a  vegeta- 
tive index)  and  the  crop  moisture  index  (CMI ) in  the  U.S.  Great  Plains 
(Exhibit  1,E).  The  Green  Number  was  used  successfully  in  the  Large 
Area  Crop  Inventory  Experiment  (LACIE)  to  monitor  and  assess  drought 
conditions  in  the  U.S.  and  foreign  areas. 

An  earlier  study  by  the  OCAD  found  a high  correlation  between  a number 
of  Vi's  computed  over  spring  wheat  fields  indicating  the  Vi's  measure 
greenness  similarly  (Exhibit  1,F). 

1 . 3 PURPOSE 

The  purpose  of  this  study  was  to  assess  the  relationship  of  the  Landsat 
data,  in  terms  of  vegetative  indices  (Vi's)  to  plant  parameters,  such 
as  plant  height,  plant  density  and  yield. 

The  Vi's  evaluated  in  this  study  included  the: 


1) 

Ashbum  Vegetative  Index 

(AVI) 

2) 

Difference  Vegetative  Index 

Green  Veget at ive  Indexrm 

Kauth  Vegetat  ive  Index‘d 

(DVI) 

3) 

(GVI) 

4) 

(KVI) 

5) 

Leaf  Area  Index 

(DAI) 

6) 

Perpendicular  Vegetative  Index  (Band  6) 

(PVI6) 

7) 

Perpendicular  Vegetative  Index  (Band  7) 

(PVI7) 

8) 

Transformed  Vegetat ive  Index  (Band  6) 

(TVI6) 

9) 

Transformed  Vegetative  Index  (Band  7) 

(TVI 7 ) 

The  VI  transformat  ions  are  shown  in  Appendix  A. 

The  objectives  of  this  study  were: 

1.  To  identify  the  VI (s)  which  provides  the  mast  significant 
relationship  to  plaint  height,  plant  density  and  yield. 

2.  To  identify  the  time(s)  in  the  growing  season  during  which 
the  relationship  between  the  VI (s)  and  plant  height,  plant 
density,  soil  moisture  and  yield  are  the  strongest. 

3.  To  identify  the  functional  relationship  of  the  Vi's,  soil 
moisture,  and  crop  calendar  to  plant  height,  plant  density, 
soil  moisture  and  yield. 


1.4  DATA  SET 

1.4.1  Image  Data. 


IJ 


The  GVI  is  equivalent  to  the  greenness  ccmponent  of  the  Kauth  transform 
and  the  KVI  is  equivalent  to  the  Green  Number  used  during  the  LACIE. 
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1.4. 1.1  Segment  Location.  The  image  data  were  extracted  from  22  blind 
sites 2/  identified  by  the  LACIE.  Seven  of  the  segments  were 
located  in  Montana,  and  the  remaining  15  were  located  in  North 
Dakota  (Figure  1-1).  The  segments  were  stratified  into  Agro- 
physical Units  (APU)  developed  by  the  LACIE.  An  APU  is  a 
geographic  area  having  similar  agricultural  and  physical 
characteristics  such  as  similar  climate,  soils,  opography  and 
agricultural  density. 

1.4. 1.2  Observations . Located  within  each  LACIE  blind  site  were  15 
spring  and/or  winter  wheat  fields  for  which  ground  observed 
information  such  as  plant  height,  plant  density  and  yield  were 
col lected. 

The  fields  in  this  study  varied  in  size  frcm  9 to  405  acres, 
averaging  41  acres.  Although  15  ground  truth  fields  were  avail- 
able, small  fields,  clouds,  haze  and  mi sregi strat ion  of  the 
image  data  limited  the  nurber  available  for  analysis  to  approx- 
imately 10  fields  per  blind  site. 

Landsat  data  was  acquired  between  April  and  September  of  1978 
and  yielded  103  useable  acquisitions.  Most  of  the  acquisitions 
were  collected  during  the  early  and  late  parts  of  the  growing 
season  as  shown  below. 


Grcwth  Stage 

% of  Obse 

Planting 

9 

Ti 1 lering 

21 

Stem  Extension 

3 

Heading 

5 

Flowering 

10 

Ripening 

30 

Harvest 

22 

The  VI  statistical  mean  was  ccmputed  for  each  field  within  each 
blind  site.  A total  of  980  observations  were  defined  frcm  the 
103  useable  acquisitions. 

1.4. 1.3  Image  Corrections.  Preprocessing  of  the  digital  data  was  exer- 
cised prior  to  the  calculation  of  the  Vi's.  Both  Landsat  2 
and  3 data  were  used  in  the  study.  Landsat  3 data  were 
calibrated  to  the  Landsat  2 data.  Channels  1-4  were  multi- 
plied by  the  calibration  values  of  1.161,  1.230,  1.246  and 
1.062,  respectively.  These  values  were  developed  by  the 
Lockheed  Electronics  Ccmpany  and  are  consistent  with  those 
suggested  by  the  Environmental  Research  Institute  of  Michican 
(ERIM)  in  a recent  study  "Landsat-3  to  Landsat-2  Calibration 
Transformation",  (NASA  Contract  NAS  9-15476). 


2 J 

A LACIE  blind  site  is  a 5x6  nautical  mile  sample  segment  for  which 
ground  observed  data  such  as  field  identifications  were  collected. 
Blind  site  information  was  used  by  the  LACIE  to  assess  the  accuracy 
of  crop  area  classification  results  achieved  by  the  LACIE  analysts. 
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FIGURE  1-1  AGROPHYSICAL  UNIT  AND  SEGMENT  LOCATIONS 
IN  MONTANA  AND  NORTH  DAKOTA 
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The  image  data  were  also  sun  angle  corrected  to  a 39°  solar 
zenith  angle.  Sun  angle  correction  compensates  for  differ- 
ences in  digital  radiance  due  to  imagery  acquired  at  different 
latitudes  and  at  different  t imes  during  the  growing  season. 

Prior  to  haze  correction,  the  data  were  autcmat ical ly  screened 
for  bad  data.  The  screening  algorithm  flags  and  removes  from 
the  haze  correction  algorithm  those  data  elements  (pixels)  that 
are  garbled,  clouds,  water  and  cloud  shadows.  After  removing 
"bad"  data,  the  haze  correction  algori thm  more  accurately  compen- 
sates for  the  effects  of  atmospheric  attenuation  due  to  haze. 

The  screening  and  haze  algorithms  used  were  the  XSTAR  developed 
by  ERIM.  The  XSTAR  correction  parameters  were  developed  for 
Landsat  2 data  and  adjusted  to  a solar  zenith  angle  of  39°.  The 
ERIM  study,  previously  cited,  also  indicated  that  the  screening 
and  haze  correction  functions  developed  for  Landsat  2 data 
operated  reasonably  well  with  calibrated  Landsat  3 data. 

1.4.2  Ground  Observed  Data.  Ground  observed  data  for  the  15  wheat  fields 
consisted  of  wheat  type,  plant  height  (inches),  ground  cover 
(percentage),  growth  stage  (Feeke's  scale),  yield  (bushels  per 
acre),  and  field  carments.  During  the  growing  season  experienced 
county  Agricultural  Stabilization  and  Conservation  Service  (ASCS) 
personnel  visited  the  blind  sites  on  satellite  overpass  days  for 
ground  data  collection.  Ground  data  information  recorded  is  self- 
explanatory  except  for  the  ground  cover  and  growth  stage  data. 

Ground  cover,  or  plant  density  was  recorded  as  codes  1-5.  Each 
code  represented  a 20  percent  range  in  ground  cover.  Growth  stage 
data  were  recorded  according  to  the  stages  of  wheat  development 
(Figure  1-2)  devised  by  W.  Feekes  of  the  Netherlands.  The  accuracy 
of  the  collected  ground  data  can  not  be  ascertained,  although 
highly  accurate  information  would  have  required  a good  agronomist. 
Study  of  the  yield  data  suggest  that  it  was  generally  obtained  by 
visual  observation,  however,  farmer  estimates  were  available  for 
seme  fields. 

1.4.3  Soil  Moisture  Data.  Meteorological  data  collected  at  weather 
stations  ident i f i ed  near  each  blind  site  was  used  to  estimate 
soil  moisture.  Temperature  and  precipitation  data  were  obtained 
from  April  through  September  of  1978  and  provided  the  input  to 
the  Two  Layer  Soil  Moisture  Model  used  by  the  National  Oceano- 
graphic and  Atmospheric  Administration  (NOAA)  to  compute  their 
Crop  Moisture  Index  (CMI).  The  model  uses  long  term  temperatures, 
current  temperatures  and  precipitation,  and  potential  water  holding 
capacity  of  the  soils  to  estimate  soil  moisture  levels.  The  surface 
and  subsurface  soil  moisture  estimates  obtained  from  the  model  are 
not  field  specific,  but  are  estimates  for  the  weather  stations 
located  in  close  proximity  to  the  blind  sites.  Field  level  soil 
moisture  values  would  have  been  preferred  for  this  study  since 

the  spectral  and  ground  observed  data  are  both  at  field  level. 
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FIGURE  1-2  FEEKE 'S  GROWTH  STAGE  SCALE  FOR  WHEAT 
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notes) 
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(wheat) 


Stage 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

10.1 


11.1 

11.2 

11.3 

11.4 


One  shoot  (number  of  leaves  can  be  added)  = " brairding  " 

Beginning  of  tillering 

Tillers  formed,  leaves  often  twisted  spirally.  In  some  varieties  of 
winter  wheats,  plants  may  be  “ creeping  ” or  prostrate  V 

Beginning  of  the  erection  of  the  pseudo-stem,  leaf  sheaths  beginning 
to  lengthen 

Pseudo-stem  (formed  by  sheaths  oflcavcs)  strongly  erected 
First  node  of  stem  visible  at  base  of  shoot 
Second  node  of  stem  formed,  next-to-last  leaf  just  visible 
Last  leaf  visible,  but  still  rolled  up,  car  beginning  to  swell 
Ligulc  of  last  leaf  just  visible 

Sheath  of  last  leaf  completely  grown  out,  ear  swollen  but  not  yet 
visible 

First  cars  just  visible  (awns  just  showing  in  barley,  ear  escaping 
through  split  of  shcatfi  in  wheat  or  oats) 

10.2  Ouartcr  of  heading  process  completed 

10.3  Half  ofheading  process  completed 

10.4  Three-quarters  of  heading  process  completed 

10.5  All  cars  out  of  sheath 

10.5.1  Beginning  of  flowering  (wheat) 

10.5.2  Flowering  complete  to  top  of  ear 

10.5.3  Flowering  over  at  base  of  car 

10.5.4  Flowering  over,  kernel  Watery  ripe 
Milky  ripe 

Mealy  ripe,  contents  of  kernel  soft  but  dry 
Kernel  hard  (difficult  to  divide  by  thumb-nail) 

Ripe  for  cutting.  Straw  dead 


Tn  if  juno 


Stem 

Extension 


Headino 


Feowejuno 

(Wheat) 
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1.5  APPROACH 

The  nominal  relationship  between  a VI  and  growth  stage  is  shown  as 
Figure  1-3.  The  tenporal  profile  rises  from  time  of  emergence,  peaks 
at  or  around  heading  and  then  drops  off  sharply  to  harvest.  The  shape 
of  the  profile  for  each  VI  is  similar.  Those  factors  affecting  the 
height  and/or  slope  of  the  temporal  profile  other  than  the  growth  stage 
factor  are  variables  such  as  crop  condition,  varietal  differences  and 
cultural  practices.  Assuming  the  growth  stage  information  obtained 
from  the  ground  truth  data  is  fairly  accurate,  the  analysis  can  then 
proceed  to  identify  those  factors  affecting  the  height  and/or  slope  of 
the  prof i le. 

The  relationship  between  a VI  and  plant  density  is  positive  (i.e.,  as 
the  biomass  increases  in  density,  the  value  of  the  VI  also  increases). 
Nominally,  the  denser  the  biomass  at  a given  growth  stage,  the  better 
the  condition  of  the  crop  and  its  yield.  There  are  exceptions  to  this 
relationship,  but  generally  the  preceding  statement  is  true.  One  such 
exception  is  the  case  when  insects  attack  the  fruit  of  the  crop,  but 
do  not  affect  the  vegetative  growth.  The  density  of  the  green  bicmass 
remains  unaffected  and  is  therefore,  undetectable  by  Landsat. 

For  purposes  of  the  analysis,  the  tenporal  profile  was  subdivided  into 
seven  distinct  intervals  as  shown  below.  The  tillering  through  ripen- 
ing growth  stage  intervals  are  based  on  the  Feekes  scale,  while  the 
planting  and  harvest  intervals  were  added  for  purposes  of  this  analysis 
and  were  given  the  0.0  and  12.0  growth  scale  definitions,  respectively. 


Interval 

Growth  St 

ages 

Label 

1 

0.0 

Plant ing 

2 

1.0 

5.0 

Til lering 

3 

6.0 

- 10.0 

Stem  Extension 

4 

10.1 

- 10.5 

Heading 

5 

10.51 

- 10.54 

Flowering 

6 

11.1 

- 11.4 

Ripening 

7 

12.0 

Harvest 

Rather  than  define  the  profile  by  a single  function,  the  profile  was 
defined  as  seven  individual  ranges  or  functions.  Several  transfor- 
mations were  tested  to  depict  the  planting  through  ripening /harvest 
relationship  between  a VI  and  growth  stage,  but  they  were  not  consistent. 
It  was  hypothesized  the  relationship  between  a VI  and  crop  condition 
parameters  (such  as  yield)  within  a growth  stage  interval  would  approx- 
imate a linear  relationship.  Changes  in  the  slope  of  the  profile 
within  a given  interval  are  minimal  compared  to  the  changes  in  the  slope 
over  the  entire  probile.  The  analysis  identified  the  relationships 
between  the  Vi's  and  crop  condition  critireria  for  each  of  the  seven 
intervals.  A major  constraint  of  the  analysis  was  the  uneven  dis- 
tribution of  the  data  within  the  seven  intervals.  Those  intervals 
within  which  the  strongest  relationships  were  hypothesized  were  the 
stem  extension,  heading  and  flowering  intervals.  Unfortunately,  they 
accounted  for  only  18  percent  of  the  data.  The  uneven  distribution 
of  the  data  across  the  growth  stage  intervals  was  a severe  limitation 
in  this  study. 
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FIGURE  1-3  RELATIONSHIP  BETWEEN  GROWTH  STAGE  AND  A VI 
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Part  2.0  reports  the  results  of  the  evaluation  of  general  relation- 
ships between  the  Vi's  and  ground  observed  data  at  each  of  the  seven 
intervals.  The  surface  and  subsurface  soil  moisture  levels  were  added 
to  the  data  set  to  iiqjrove  the  likelihood  of  identifying  a strong  and 
repeatable  relationship  between  yield  and  the  indices.  It  was  felt  the 
Vi's,  along  with  surface  and  subsurface  soil  moisture  parameters,  would 
account  for  a larger  percentage  of  the  variation  in  yield,  than  would 
be  achieved  by  using  the  indices  alone. 

Correlation  analysis  was  used  to  determine  the  strength  of  the  relation- 
ships between  the  Vi's  and  crop  data  at  each  growth  stage  interval. 
Specifically,  all  the  VI ' s were  correlated  to  each  other,  plant  height, 
plant  density,  growth  stage,  yield,  and  surface  and  subsurface  soil 
moisture.  A matrix  of  correlation  coefficients  and  significance  levels 
at  each  growth  stage  interval  is  presented  in  Part  2.0.  It  was 
hypothesized  that  the  relationship  between  the  VI ' s and  plant  density, 
plant  height  and  yield  would  strengthen  as  the  plant  progressed  through 
its  growth  cycle  reaching  the  strongest  relationship  at  heading  and 
flowering,  and  weakening  at  ripening  and  harvest  time. 

Additionally,  analysis  of  the  coefficient  of  variation  (CV)  of  a VI 
at  a given  growth  stage  was  implemented  to  determine  if  the  CV  of  a 
VI  is  reduced  by  stratifying  the  observations  into  predefined  yield 
intervals  versus  the  CV  of  a VI  across  all  yields  (unstratified).  The 
observat ions  were  stratified  into  five  bushel /acre  yield  intervals 
from  10-14  bu/ac  to  60-64  bu/ac.  Yields  ranged  frcm  a low  of  12  bu/ac 
to  a high  of  63  bu/ac  and  averaged  34.5  bu/ac. 

The  data  were  further  analyzed  to  determine  the  effects  of  field  size, 
wheat  type  (winter  or  spring)  and  APU  location  on  the  relationships 
between  the  Vi's  and  the  crop  data.  Correlation  coefficients  were 
ccmputed  for  t^ch  of  the  three  field  size  ranges:  greater  than  30,  40, 
and  50  pixels  . It  was  felt  the  relationship  would  strengthen  as  the 
minimum  field  size  was  increased.  Increasing  the  rninirmm  field  size 
would  reduce  the  effect  of  edge  or  border  pixels  on  the  VI  statistical 
mean  computed  for  each  field.  Correlation  coefficients  were  computed 
separately  for  those  winter  wheat  and  spring  wheat  observations.  This 
analysis  would  indicate  the  effects  of  wheat  type  on  the  relationships 
of  the  Vi's  to  plant  density,  plant  height,  and  yield.  Correlation 
coefficients  wore  also  computed  for  those  observations  contained  in 
each  of  four  APU's.  It  was  hypothesized  that  relationships  would 
strengthen  as  effects  due  to  soils,  topography  and  climate  were  mini- 
mized. An  APU  is  a geographic  area  having  similar  physical  and 
agricultural  characteristics,  such  as  similar  soils,  topography, 
agricultural  density  and  climate. 

Part  3.0  presents  the  results  of  the  regression  analysis  performed  on 
the  data  set.  A stepwdse  regression  procedure  was  performed  at  each 
of  the  seven  growth  stage  intervals  for  each  of  nine  Vi's  using  yield, 
plant  density  and  plant  height  as  the  dependent  variables  and  the 
Vi's,  surface  and  subsurface  soil  moisture,  growth  stage  and  a number 


A pixel  or  picture  element  is  the  smallest  area  for 
measurements  are  collected  by  Landsat.  A Landsat  pixel 
1.1  acres  in  size. 
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of  interactive  terms  as  the  independent  variables.  Relationships  were 
carpared  for  each  VI  by  growth  stage  interval  in  an  effort  to  define 
the  time  in  the  growing  season  when  the  strongest  relationship  occurs. 
Additionally,  the  regression  results  for  all  the  Vi's  were  compared  to 
one  another  by  growth  stage  interval  in  an  effort  to  identify  the  VI 
that  provided  the  strongest  relationship  to  yield.  The  coefficients 
of  determination  (R  ),  standard  deviations  and  other  associated 
statistics  were  computed  and  are  presented  in  Part  3.0.  The  standard 
deviation  (i.e.,  square  root  of  the  mean  square  error  (N/  M5E)  of  the 
difference  between  the  observed  and  estimated  yield  is  the  principal 
statistic  used  in  the  comparison  of  the  yield  relationships.  Each 
regression  relationship  determined  at  a given  growth  stage  interval 
for  a VI  is  based  on  the  identical  data  set.  The  only  element  \diich 
varied  was  the  particular  VI  entered  into  the  analysis.  Similar 
interactive  terras  were  defined  for  each  VI  and  entered  into  the  step- 
wise regression  procedure.  An  interactive  term  is  the  cross-product 
between  two  or  more  variables,  such  as  the  cross-product  of  a VI  and 
soil  moisture  or  growth  stage  (e.g.,  AVI*Soil  Moisture). 

A 10  percent  level  of  significance  was  introduced  as  a constraint 
for  the  entry  of  a variable  into  the  regression  relationship. 

The  Statistical  Analysis  System  (SAS)  software  routines  were  used 
to  analyze  the  data.  SAS  is  an  integrated  system  for  data  manage- 
ment and  statistical  analyses. 

1.6  RECOMMENDATIONS  FDR  FIXTURE  RESEARCH 

° The  acquisition  of  a complete  and  accurate  set  of  ground  observed 
and  Lands at  data  to  support  both  static  and  dynamic  analyses  of  the 
relationships  between  crop  condition  parameters  such  as  biomass 
and  yield,  and  the  Landsat  data  is  reccrrmended. 

° An  in-depth  analysis  of  vdiat  effects  soils,  geographic 
location  and/or  other  factors  have  on  these  relationships 
is  strongly  reccmrended. 

° Actual  field  level  measurements  of  soil  moisture  is  reccmmended 
rather  than  using  estimates  obtained  frcm  a soil  moisture 
model.  There  is  a strong  indication  that  the  vegetative 
indices,  crop  calendar  and  soil  moisture  can  be  used  effectively 
to  measure  crop  condition. 
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FIGURE  2-1  LAI  VS.  GROWTH  STAGE 
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intervals  ranged  frcm  10-14  bu/ac  to  60-64  bu/ac.  The  average 
CV  for  the  eleven  yield  intervals  identified  ranged  frcm  a 
high  of  400  at  emergence  to  a low  of  4.44  at  heading  and 
averaged  53.  Stratification  did  not  result  in  an  appreciable 
reduction  in  the  average  CV  for  the  eleven  yield  intervals. 
However,  during  the  stem  extension  to  flowering  intervals  the 
CV  was  reduced  frcm  28  to  22,  a decrease  of  21  percent.  The 
variation  in  the  LAI  after  stratification  by  yield  remained 
high. 

Information  on  other  factors  contributing  to  the  large  vari- 
ation in  the  LAI  such  as  cropping  practices  was  not  available 
and  prevented  any  further  conclusive  analysis. 

2.2.2  Yield  and  Vi's. 


2. 2. 2.1  Over  Entire  Data  Set.  The  approach  used  in  this  study  was 
to  identify  the  relationship  between  yield  and  a VI(s)  at 
a point  in  time  (static  analysis)  versus  an  approach  defin- 
ing the  relationship  over  time  (dynamic  analysis).  The 
lack  of  timely  data  for  a number  of  fields  across  the  crop 
calendar  prevented  a more  dynamic  analysis.  A more  dynamic 
approach  would  have  been  preferred,  because  the  yield  of  a 
plant  is  not  determined  at  any  one  time,  but  is  a combi- 
nation of  many  physiological  and  envi rormental  factors  that 
interact  during  the  entire  growth  cycle.  However,  the  static 
approach  allowed  the  analysis  to  identify  the  times  in  the 
crop  calendar  during  vhich  yield  and  a VI  were  strongly 
correlated. 

Figures  2-2  to  2-8  shew  the  relationship  between  yield  and 
the  LAI  at  each  of  seven  growth  stage  intervals.  A single 
growth  stage  within  an  interval  was  used  to  generate  the 
plot  in  an  effort  to  minimize  the  variation  in  the  LAI 
wihin  a growth  stage  interval.  For  example,  growth  stage 
10.0  of  the  stem  extension  growth  stage  interval  was  used 
to  generate  Figure  2.4  (the  stem  extension  interval  extends 
from  growth  stages  6.0  - 10.0  on  the  Feeke's  scale). 

The  relationship  between  yield  and  the  LAI  was  highly  vari- 
able during  each  growth  stage  interval.  Hie  greatest 
variations  occurred  during  the  planting,  tillering,  ripen- 
ing and  harvest  growth  stage  intervals,  vhile  the  least 
variation  occurred  during  the  stun  extension,  heading  and 
flowering  growth  stage  intervals.  The  LAI  alone  was  not 
sufficient  for  explaining  a large  percentage  of  the  variation 
in  yield. 

Correlation  tables  were  generated  for  each  of  the  seven 
growth  stage  intervals  in  an  effort  to  identify  the  strength 
of  the  interrelationship  between  the  spectral  data  and 
the  ground  observed  data.  Correlation  Tables  2-1  to  2-7 
present  these  interrelationships  (correlation  coefficients 
and  significance  levels)  of  15  variables  including  nine 
Vi's,  three  ground  observed  elements  including  plant  height, 
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FIGURE  2-2  YIELD  VS  LAI  At  PLANTING  .. 

PLOT  OF  YIELO°LAI  LEGEND:  A = 1 OPS*  8=2  0B5»  ETC 
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FIGURE  2-3  YIELD  VS, LAI  AT  TILLERING 
PLOT  OF  Y I ELD*L A I LEGEND:  A = 1 OBS»  B = 2 OBS»  ETC 
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FIGURE  2-4  YIELD  VS/ LAI  AT  STEM  EXTENSION 
PLOT  OF  Y I ELD°L A I LEGEND:  A = 1 OBSt  B = 2 OBS»  ETC 
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FIGURE  2-5  YIELD  VS,LAI  AT  HEADING 
PLOT  OF  Y I ELD«L A I LEGEND:  A = 1 OBS  * B = 2 OBS » ETC 
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FIGURE  2-6  YIELD  VS,  LAI  AT  FLOWERING 

OF  YIELO*LAI  LEGEND:  A = 1 ORS»  B = Z ORSt  ETC 
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FIGURE  2-7  YIELD  VS,  LAI  AT  RIPENING 
YlEin°LAI  LEGEND:  A = 1 ORS»  B = 2 OBS»  ETC 
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TABLE  2.1  CORRELATIONS  AND  SIGNIFICANCE  LEVELS  AT  PLANTING 
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-^YLD  = YIELD 

-^SMl  = SURFACE  SOIL  MOISTURE 
-^SM2  = SUBSURFACE  SOIL  MOISTURE 
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TABLE  2.2  CORRELATIONS  AND  SIGNIFICANCE  LEVELS  AT  TILLERING 
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T — 1 

1 — 1 
o 
o 
o 

o 

i — ! 

1 

r^ 

uo 

r— 1 

'T 

o 

CO 

o 

aO 

m 

O' 

o 

1 

m 

as 

os 

1 — 1 

PVI7 

o 

o 

1 — I 

aO 

1 — I 

o 

o 

o 

o 

m 

1 — 1 
o 
o 
o 

v£> 

m 

1 — 1 
o 
o 
o 

CN 

1 

i — i 

o 

o 

o 

m 

o 

i 
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OA 

O' 

CN 

1 

T— 1 

o 

o 

o 

M 

> 

o 

o 

1 — 1 

D 

OA 

i — i 

o 

o 

o 

o 

I — 1 

o 

o 

o 

1 — 1 
IT) 

1 — 1 
o 
o 
o 

m 

1 — I 

o 

o 

o 

m 

Cs] 

r 

cn 

o 

o 

o 

o 

1 

t-H 

o 

CO 

CO 

co 

CO 

I- 

1 — 1 
o 
o 
o 

Pj 

o 

o 

1 — 1 

OA 

CO 

i — i 

o 

o 

o 

r^- 

OA 

1 — 1 
o 
o 
o 

co 

vD 

1 — I 

o 

o 

o 

AO 

1 — 1 

o 

o 

o 

vO 

1 — 1 
o 
o 
o 

m 

CN 

cn 

o 

o 

o 

o 

NS 

in 

00 

m 

CN 

1 

CN 

o 

o 

o 

M 

> 

H 

TVI6 

O 

o 

OA 

i — 1 

O 

o 

o 

CO 

OA 

1 — 1 
o 
o 
o 

OA 

1 — I 

o 

o 

o 

AO 

I — 1 

o 

o 

o 

00 

*<T 

1 — 1 
o 
o 
o 

m 

r— 1 

o 

o 

o 

o 

CO 

1 — 1 
o 
o 
o 

co 

o 

1 — 1 
m 

CN 

AD 

CO 

CO 

1 — 1 
o 
o 
o 

r- 

H 

M 

> 

OA 

i — 1 

O 

o 

o 

OA 

I — 1 

o 

o 

o 

AO 

OA 

1 — 1 
o 
o 
o 

OO 

OA 

1 — 1 
o 
o 
o 

c^ 

N3 

1 — i 

o 

o 

o 

t — 1 

m 

1 — I 

o 

o 

o 

NO 

m 

1 — 1 
o 
o 
o 

vO 

CN 

1 

CN 

o 

o 

o 

1 — 
o 

1 

m 

o 

m 

cn 

r^. 

CN 

1 

1 — i 

o 

c 

<3 

o 

o 

• 

i-H 

u 

o 

o 

1 — 1 

— 1 

o 

o 

o 

OA 

1 — 1 
o 
o 
o 

r^. 

OA 

1 — 1 
o 
o 
o 

N3 

O' 

! 1 

o 

o 

o 

00 

OA 

1 — 1 
o 
o 
o 

O' 

' D 

1 — 1 
o 
o 
o 

i — i 
m 

1 — 1 
o 
o 
o 

NS 

m 

1 — 1 
o 
o 
o 

aO 

CN 

CN 

o 

o 

o 

o 

in 

o 

m 

cn 

- 

CN 

iH 

O 

.c 

o 

M 

> 

o 

o 

o 

r— 

o 

o 

1 — 1 

aD 

OA 

— 1 
o 
o 
o 

NO 

CT-i 

— 1 
o 
o 
o 

O' 

1 — { 
o 
o 
o 

OA 

1 — 1 
o 
o 
o 

1 1 

OA 

1 — I 

o 

o 

o 

m 

as 

I — 1 

o 

o 

o 

<r 

aO 

1 — 1 
o 
o 
o 

m 

1 — 1 
o 
o 
o 

00 

1 — 1 
o 
o 

3 

o 

CO 

l — 1 

o 

o 

o 

CO 

O 

-d- 

CN 

00 

AD 

CN 

CO 

t— ) 

c 

c 

o 

LAI 

C 

3 

m 

OA 

— I 

o 

o 

o 

oa 

OA 

— 1 

3 

o 

o 

OA 

OA 

r“l 

3 

O 

O 

OA 

1 — I 

o 

o 

o 

OA 

1 — I 

o 

o 

o 

m 

O' 

1 — 1 
o 
o 
o 

OA 

OA 

1 — 1 
o 
o 
o 

N3 

ns 

1 — i 

o 

o 

o 

oa 

'Cf 

.0001 

m 

m 

t — 1 

o 

o 

o 

aO 

CN 

CN 

o 

o 

o 

m 

o 

OA 

m 

OA 

O'. 

CN 

I — ] 

o 

o 

o 

i — i 
> 

o 

r— 1 

Q 

AVI 

o 

o 

• 

t — 1 

30 

i — 1 

O 

O 

o 

<T 

OA 

— i 

o 

o 

o 

00 

OA 

— 1 
3 
o 

3 

oo 

— ! 
O 
3 
O 

<r 

OA 

T — 1 

o 

o 

o 

00 

OA 

1 — I 

o 

o 

o 

OA 

00 

— 1 
o 
o 
o 

NO 

O' 

i — ! 

o 

o 

o 

CO 

aO 

I — 1 

o 

o 

o 

r-. 

1 — 1 
o 
o 
o 

1 — i 

o 

o 

o 

CN 

CN 

1 — 1 
CN 
O 

o 

o 

OA 

OA 

CN 

CO 

OJ 

1 — 1 

1 — 1 
o 
o 

M 

> 

< 

M 

> 

Q 

1— 1 

M 

> 

o 

M 

> 

AD 

M 

£ 

M 

> 

H 

PVI6 

i — i 
> 

PH 

PLTH 

PLTD 

GROW 

YIELD 

f 1 

C/0 

CN 

S 

CO 

SEP  1 8 1979 
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'LTH  = PLANT  HEIGHT 
LTD  = PLANT  DENSITY 
■ROW  = GROWTH  STAGE 
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TABLE  2^3  CORRELATIONS  AND  SIGNIFICANCE  LEVELS  AT  STEM  EXTENSION 


SM2 

f 

I.  00 

4— 

SMI 

o 

o 

rH 

CN 

3 

1 

rH 

<r 

f\ 

o 

r 

YLD 

ro 

CO 

i — 

CO 

3 

3 

NT 

m 

i — i 

00 

i — 1 

3 

H 

i 

GROW 

O 

3 

tH 

CN 

CM 

yD 

CO 

CM 

CM 

i — i 
CM 

<r 

CO 

CM 

co 

co 

o 

o 

PLTD 

— 

m 

r\ 

T 

o 

3 

3 

O 

3 

3 

CN 

3 

m 

H 

3 

3 

3 

3 

C 

C 

D 

3 

H 

o 

CM 

I 

r 

PLTH 

O 

o 

t — 1 

1 — l 

m 

3 

ro 

3 

3 

CM 

n- 

i — 1 

3 

3 

3 

r- 

CO 

05 

CO 

o 

3 

3 

YD 

CM 

YD 

3 

3 

H 

CO 

<r 

yD 

PVI7 

O 

O 

i — 1 

<r 

00 

i — ! 

O 

3 

O 

3 

cO 

i — I 
O 
3 
3 

05 

yD 

1 — l 
3 
3 
O 

yD 

CO 

00 

3 

o 

i — 1 

O 

1 

GY 

MT 

Or 

co 

i — 1 

1 

.4620 

PVI6 

O 

0 

1  1 

3 

O'! 

i — 1 

3 

3 

3 

CM 

00 

.0001 

yD 

t — I 
O 

o 

o 

05 

yD 

.0001  1 

o 

CO 

CO 

3 

3 

o 

o 

rH 

OY 

i — 1 
00 

o 

3 

00 

r-* 

CM 

TVI7 

O 

3 

i — 1 

3 

3 

i — i 

O 

o 

3 

00 

O'. 

i — 1 

O 

o 

3 

CO 

rH 

o 

3 

3 

00 

3 

i — ! 

3 

O 

O 

i — 1 

n- 

1 — 1 
3 
3 
O 

o 

.0248 

<■ 

o 

3 

3 

3 

00 

3 
t — l 

1 

.3841 

vO 

CO 

CO 

i — ! 

3 

3 

3 

3 

3 

1 1 

O 

3 

3 

n- 

On 

rH 

3 

3 

3 

0-* 

i— 1 

O 

o 

o 

1 1 

T 1 

o 

3 

3 

CO 

i — 1 

3 

O 

3 

MO 

CO 

o 

3 

<r 

3 

CM 

O 

3 

3 

3 

00 

3 
< — 1 

1 

3 

CO 

3 

3 

H 

£ 

o 

o 

1 — 1 

H 

M 

o 

o 

i — i 

00 

3 

T 1 

o 

o 

o 

00 
a '» 

t — 1 

3 

3 

3 

3 

3 

i — i 

3 

3 

3 

3 

o •» 

1000* 

CM 

00 

.0001 

00 

3 

rH 

3 

3 

3 

O 

rH 

3 

3 

O 

3 

CO 

i — i 

Ln 

o 

CN 

3 

CO 

O 

3 

00 

1 — 1 

1 

YD 

05 

CO 

CO 

M 

> 

o 

o 

I — 1 

o 

o 

1 — 1 

i — 1 

o 

o 

o 

00 

3 

1 — I 

o 

o 

o 

00 

3 

i — 1 

3 

3 

O 

(3 

3 

i — 1 

3 

3 

3 

3 

3 

1 — 1 
3 

o 

3 

3 

OO 

i — i 

3 

O 

3 

CO 

\D 

i — 1 

3 

O 

3 

o 

i — i 

3 

O 

3 

3 

CO 

■ — 1 

m 

o 

CN 

O 

00 

~cr 

3 

3 

00 

rH 

yD 

05 

CO 

CO 

o 

M 

<j 

o 

o 

rH 

00 

3 

.0001 

oo 

3 

.0001 

3 

3 

— 1 
o 
o 
o 

CO 

CO 

.0001 

vO 

3 

i — 1 

3 

3 

3 

CO 

i — 1 

3 

3 

3 

yD 

r->* 

i — 1 

O 

o 

o 

CO 

3 

i — 1 

3 

3 

3 

rH 

i — ! 

3 

3 

3 

00 

ro 

<r 

co 

o 

3 

O 

00 

3 

3 

3 

CO 

CM 

f 

05 

o 

CM 

CM 

*— 1 

DVI 

i 

o 

o 

i — ! 

CO 

C5 

— 1 
o 
o 
o 

Or 

05 

— 1 
o 
o 
o 

3 

3 

— i 

o 

o 

o 

3 

O'! 

— 1 
3 

o 

o 

CO 

I — 1 

o 

3 

O 

00 

3 

— 1 
3 
O 
O 

o 

3 

i — 1 

O 

3 

3 

*<r 

00 

I — 1 

o 

3 

O 

Ln 

yd 

— i 

3 

3 

O 

oo 

3 

1 1 

3 

3 

3 

yD 

CO 

3 

3 

<r 

3 

I { 

o 

<r 

i — i 

05 

CO 

rH 

H 

O 

AVI 

o 

o 

rH 

o 

o 

1 — { 

i — i 

o 

o 

o 

00 

3 

1000* 

CO 

CO 

1000* 

CO 

CO 

— 1 
o 
o 
o 

DO 

CO 

— 1 

3 

3 

3 

CO 

CO 

t — 1 

3 

3 

3 

00 

CO 

r— 1 

3 

3 

3 

3 

3 

— i 

O 

3 

3 

co 

oo 

i — 1 

3 

O 

3 

m 

■JD 

— 1 
3 
3 
3 

3 

3 

.0001 

CO 

.0385 

3 

3 

YD 

05 

rH 

1 

Ln 

M 

< 

1 — 1 
> 

Q 

M 

< 

3 

M 

> 

o 

M 

> 

TVI6 

o- 

M 

£ 

H 

> 

P-. 

3 

M 

> 

P-l 

E- 

K 

CH 

GROW 

YIELD 

r 1 

s 

CO 

3 

s 

CD 

SEP  1 8 1379 
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TABLE  2-4  CORRELATIONS  AND  SIGNIFICANCE  LEVELS  AT  HEADING 


SM2 

j 

1 

f 

i 

o 

o 

— 1 

i 

I 

SMI 

- 

o 

o 

1 — 1 

c 

c 

.9905 

r— 

1 

YLD 

— 

LT 

c 

CN 

co 

O' 

LT 

i — i 

o 

o 

o 

- 

o 

— 1 

1 

GROW 

o 

o 

T— 1 

r^« 

On 

1 

CN 

a 

c- 

c 

cr 

c 

vO 

CO 

T— 

m 

O- 

On] 

1 

so 

00 

co 

o 

t 

PVI7 

o 

o 

1 — 1 

o 
1 — 

1 

G> 

CO 

<r 

CO 

in 

T — 

CD 

O 

o 

LO 

CN 

o 

o 

nO 

o 

Os] 

rH 

O- 

o 

PVI6 

L_ 

o 

o 

I — 1 

00 

CT5 

i — 

o 

o 

o 

as 

rH 

1 

o 

vo 

r— 1 

G\ 
l n 

1 — 1 
o 
o 
o 

in 

OS] 

.0660 

co 

i — l 
i — l 

o 

TVI7 

o 

o 

1 — 1 

<r 

OA 

i — 1 

O 

o 

o 

o** 

ON 

1 ! 

o 

o 

o 

on 

o 

LO 

m 

o 

m 

o 

NO 

rH 1 

o 

o 

o 

i — 1 

CN] 

o 

NT 
Osl 
i — 1 

05 

CM 

i .0341 

r 

I 

TVI6 

o 

o 

On 

on 

i — i 

o 

o 

o 

m 

on 

1 — 1 
o 
o 
o 

UO 

cn 

1 — 1 
c 

O' 

o 

CO 

i — 1 

CO 

On] 

CO 

ro 

SO 

1 — 1 
o 
o 
o 

o 

CM 

Os] 

r— 1 

CN 

CO 

rH 

O- 

r~H 

O 

T— 

-i 

H 

o 

o 

00 

Os 

i — i 

o 

o 

o 

1 — 1 
o 
o 
o 

On 

on 

.0001 

G^ 

o\ 

i — ! 

o 

o 

o 

LT) 
i — 1 

O 

r^. 

C\] 

05 

m 

1 — 1 
o 
o 
o 

Os] 

m 

in 

o 

i — 1 

CO 

00 

on 

i — i 

o 

PS 

r— 

H 

o 

o 

i — i 

i — 1 

O 

o 

o 

CO 

as 

1 — 1 
o 
o 
o 

CT\ 

1 — 1 
o 
o 
o 

as 

as 

i — i 

o 

o 

o 

a\ 

1 — i 

o 

o 

o 

LO 
r— 1 

o 

On] 

05 

m 

1 — 1 
o 
o 
o 

CN 

uo 

in 

o 

r— 1 
CO 

00 

as 

i — ! 

O 

M 

> 

; 

o 

o 

1 — 1 

o 

o 

o 

1 — 1 

m 
a s 

i — i 

o 

o 

o 

m 

on 

1 — 1 
o 
o 
o 

o\ 

1 — 1 
o 
o 
o 

00 

1 — 1 
o 
o 
o 

CO 

cn 

1 — 1 
o 
o 
o 

CO 

o> 

1 ! 

o 

o 

o 

i — 1 
i — 1 

00 

i — i 

<T 

vO 

1 — 1 
o 
o 
o 

* — { 
On] 

o 

00 

1 — 1 

rH 

On] 

CO 

ON 

m 

i — i 

o 

M 

< 

hH 

1 

o 

o 

1 — 1 

co 

C7 

i — i 

o 

o 

o 

ON 

On 

— 1 
o 
o 
o 

G\ 

— 1 
o 
o 
o 

in 

as 

1 — 1 
o 
o 
o 

ON 

— I 

o 

o 

o 

00 

as 

1 — 1 
o 
o 
o 

o 

o 

1 — 1 
o 
o 
o 

o 

1 — 1 

O' 

vO 

*<r 

co 

uo 

1 — 1 
o 
o 
o 

in 

CN] 

ON 

CTn 

m 

o 

On] 

On] 

o 

> 

Q . 

1— 1 

o 

o 

1 — I 

o 

o 

— 1 

1 — 1 
o 
o 
o 

m 

O'' 

.0001 

CTN 

CTn 

— 1 
o 
o 
o 

as 

as 

— 1 
o 
o 
o 

g\ 

— 1 
o 
o 
o 

00 

G\ 

— 1 
o 
o 
o 

CT\ 

— 1 
o 
o 
o 

o 

o 

— 1 
o 
o 
o 

05 

O 

o 

o 

LO 

<r 

in 

— I 

o 

o 

o 

m 

00 

m 

vD 

o 

N 

o 

m 

<r 

o 

1— 1 
> 

< 

H 

> 

o 

i— i 

3 

i— i 
> 
o 

M 

TV  1 6 

M 

£ 

no 

1 — 1 
> 

PL, 

M 

> 

GROW 

YIELD 

t 1 

X 

CO 

On] 

tn 

SEP  1379 
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TABLE  2-5  CORRELATIONS  AND  SIGNIFICANCE  LEVELS  AT  FLOWERING 


I 


SM2 

1 

1 

1 

I 

o 

o 

1 

i 

SMI 

[ 

- 

3) 

X 

cM 
« — i 

<r 

o 

X 

X 

r 

r 

j 

1 

YLD 

! 

o 

o 

1 — 1 

X 

X 

r 

X 

O'! 

o 

o 

00 

O' 

1 — 1 
o 
o 
o 

j 

j— 

1 

GROW 

1 

3 

3 

o- 

CM 

O' 
m 
\ — i 

o 

o 

<r 

l 

1 — 1 
o 
o 
o 

X 

o 

1 

ON 

00 

CO 

— 

1 

c 

r— 1 

PVI7 

! 

o 

o 

rH 

X 

o 

1 

r^ 

o- 

o- 

00 

X 

o 

o 

X 

o 

<3- 

o 

X 

X 

X 

o 

o 

T ( 

CM 

i — l 

o 

O' 

o 

! 

PVI6 

I 

o 

o 

1 — 1 

X 

O' 

i — ! 

o 

o 

o 

1 — 1 
o 

1 

.9098 

o 

1 — 1 

X 

X 

o 

X 

O' 

CM 

X 

X 

rH 

o 

00 

rH 

X 

X 

X 

o 

TVI7 

o 

o 

1 — 1 

X 

X 

1 — 1 
o 
o 
o 

rH 

ON 

1 — 1 
o 
o 
o 

o 

CT\ 

X 

in 

X 

NO 

CM 

X 

o 

1 — 1 
o 

O' 

CM 

X 

X 

1 — 1 
o 

1 — 1 
CO 

X 

1 — 1 
o 
o 

o 

o 

rH 

O'! 

1 — 1 
o 
o 
o 

00 

1 — 1 
o 
o 
o 

X 

X 

1000' 

o 

X 

X 

X 

X 

o 

CO 

CM 

o 

o 

i — 1 

CM 

o 

O' 

X 

o 

CM 

CO 

o 

! 1 

o 

o 

TV  If 

M 

o 

o 

rH ( 

rs! 

O' 

1 — 1 
o 
o 
o 

CO 

ON 

1 — 1 
o 
o 
o 

O' 

O' 

rH 

o 

o 

o 

On 

On 

1 — 1 
o 
o 
o 

o 

LO 

1 — 1 
CO 
00 

O' 

1 — 1 

X 

X 

X 

i — ! 

X 

X 

0121 

CM 

CM 

0260 

M 

> 

o 

o 

iH 1 

o 

o 

— 1 

1 — 1 
o 
o 
o 

CM 

ON 

1 — 1 
o 
o 
o 

CO 

ON 

.0001 

O' 

O' 

I — 1 

o 

o 

o 

O' 

O' 

1 — 1 
o 
o 
o 

X 

o 

LO 

1 — 1 

CO 

00 

o 

1 — ! 

X 

X 

X 

1 — 1 

X 

X 

0121 

CM 

CM 

o 

X 

X 

o 

LAI 

o 

o 

i — i 

X 

00 

1 — 1 
o 
o 
o 

X 

XI 

1000' 

X 

1 — 1 
o 
o 
o 

ON 

I — 1 

o 

o 

o 

1 — 1 
00 

1 — 1 
o 
o 
o 

r— 1 
00 

1 — I 

o 

o 

o 

LO 

1 — 1 

CM 

o 

co 

rH 

' — 1 
CO 

X 

1 — 1 
o 
o 

ON 

i — 1 

o 

X 

X 

o 

CO 

co 

o 

1 — 1 
o 
o 

DVI 

o 

o 

1 — 1 

— 1 

00 

i — I 

o 

o 

o 

O' 

O' 

1 — 1 
o 
o 
o 

On 

ON 

1 — 1 
o 
o 
o 

X 

X 

1 — 1 
o 
o 
o 

1 — 1 
ON 

1 — I 

o 

o 

o 

X 

O' 

1 — 1 
o 
o 
o 

o 

o 

1 — 1 
o 
o 
o 

CM 

o 

O' 

in 

X 

00 

o 

X 

X 

o 

<3- 

o 

CO 

X 

X 

o 

o 

rH 

CM 

CM 

o 

O’ 

o 

AVI 

o 

o 

i — i 

o 

o 

— 1 

.0001 

x 

00 

1 — 1 
o 
o 
o 

on 

ON 

1 — 1 
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TABLE  2-6  CORRELATIONS  AND  SIGNIFICANCE  LEVELS  AT  RIPENING 
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plant  density  and  yield,  and  two  soil  moisture  estimates 
including  surface  and  subsurface  soil  moisture  at  the 
seven  growth  stage  intervals. 

The  level  of  correlation  among  the  Vi's  was  high  throughout 
the  growth  cycle.  These  relationships  reached  a maximum  at 
stem  extension  and  heading.  The  correlation  between  Vi's 
based  on  band  6 to  those  Vi's  based  on  band  7 were  lower 
when  corrpared  to  correlations  of  the  Vi's  based  on  the  same 
near  infrared  band.  The  LAI,  TVI6,  and  PVI6  use  band  6 in 
their  calculations,  while  the  AVI,  DVI , GVI , KVI , TVI7,  and 
PVI7  use  band  7.  Generally,  the  LAI  and  TVI6  consistently 
produced  a stronger  relationship  to  yield  at  each  of  the 
growth  stage  intervals. 

The  relationship  between  the  Vi's  and  yield  and  the  other 
ground  observed  parameters  were  at  approximately  the  same 
level  of  correlation  due  to  the  high  inter-correlation 
between  the  Vi's.  The  correlations  did  vary  between  these 
variables,  but  they  were  not  significantly  different.  The 
relationship  between  yield  and  the  Vi's  was  not  expected 
to  strengthen  until  stem  extension  and  peak  until  heading. 

Tie  average  correlation  coefficients  between  the  Vi's  and 
yield  were  .36  at  stem  extension,  .58  at  heading,  and  .17 
at  flowering.  Correlation  coefficients  between  yield  and 
the  Vi's  at  planting  and  harvest  were  negative  and  averaged 
-.31  and  -.42,  respectively. 

Correlation  coefficients  between  yield  and  the  Vi's  at  tiller- 
ing were  negative  and  averaged  -.26.  The  increase  in  the  plant 
bicmass  reduced  the  correlation  coefficients  found  at  tillering. 
The  Vi's  were  no  longer  measuring  the  reflectance  from  the  soils 
primarily,  due  in  part  to  the  level  of  moisture,  but  also 
beginning  to  measure  the  reflectance  of  the  plant  canopy. 


Correlation  coefficients  between  yield  and  the  Vi's  at  ripening 
varied  from  being  slightly  negative  to  being  slightly  positive. 
The  correlation  coefficients  averaged  .01  and  indicated  a weak 
relationship  between  yield  and  the  Vi's. 

2. 2. 2. 2 Effect  of  Wheat  Type.  The  intentions  of  this  analysis  were  to 
compute  correlation  tables  similar  to  those  presented  earlier, 
independently  for  spring  wheat  and  winter  wheat  observations  at 
the  stem  extension,  heading  and  flowering  growth  stage  intervals. 
The  lack  of  a sufficient  number  of  observations  for  each  wheat 
type  at  each  growth  stage  interval  prevented  a complete  analysis 
across  the  entire  crop  calendar.  The  data  did  support  analysis 
for  both  wheat  types  at  heading,  but  not  at  stem  extension  and 
flowering. 

Stratification  by  wheat  type  did  not  affect  the  relationship 
between  yield  and  the  Vi's  at  heading  (refer  to  Appendix  B for 
Correlat ion  Tables) . The  average  correlation  coefficient  was 
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.57  for  spring  wheat  and  .61  for  winter  wheat,  compared  to  an 
average  correlation  coefficient  of  .58  for  the  entire  data  set, 
irrespective  of  wheat  type. 

Stratification  by  wheat  type  did  not  affect  the  relationship 
between  yield  and  the  Vi's  at  flowering.  The  average  correl- 
ation coefficient  at  flowering  was  .19  for  spring  wheat  compared 
to  .18  for  all  wheat.  Analysis  could  not  be  conducted  independ- 
ently for  winter  wheat  due  to  an  insufficient  nurber  of 
observations . 

Stratification  by  wheat  type  at  stem  extension  did  have  same 
effect  on  the  relationship  between  yield  and  the  Vi's.  The 
average  correlation  coefficient  was  .48  for  winter  wheat 
compared  to  .36  for  all  wheat,  an  increase  of  33  percent. 

2. 2. 2. 3 Effect  of  Field  Size.  The  purpose  of  increasing  the  minimum 
field  size  was  to  reduce  the  effects  of  border  or  edge  pixels 
on  the  VI  statistical  mean  computed  over  each  field.  The 
average  field  size  for  fields  greater  than  30  pixels  was  57 
pixels,  for  fields  greater  than  40  pixels  the  average  size  was 
78  pixels  and  for  fields  greater  than  50  pixels  the  average 
size  was  91  pixels.  The  average  size  for  all  fields  in  the 
data  set  was  41  pixels,  or  approx imat ely  45  acres. 

Generally,  the  field  size  constraint  reduced  the  strength  of 
the  relationship  between  yield  and  the  Vi's  at  stem  extension 
(refer  to  Appendix  C for  Correlation  Tables).  The  average 
correlation  coefficient  between  yield  and  the  Vi's  was  .53, 

.23,  and  .14  for  the  30,  40  and  50  pixel  field  size  constraints, 
respectively.  This  compares  wd th  the  average  correlation  co- 
efficient of  .36  for  all  the  fields  in  the  data  set.  The 
relationship  between  yield  and  the  Vi's  is  relatively  low 
and  inconsistent  during  the  stem  extension  period. 

During  heading,  the  effect  of  the  field  size  constraints  was 
to  strengthen  the  relationship  between  yield  and  the  Vi's. 

The  average  correlation  coefficient  increased  to  .73,  .72 
and  .66  for  the  30,  40  and  50  pixel  field  size  contraints, 
respectively.  This  compares  with  the  average  correlation 
coefficient  of  .58  at  heading  based  on  all  the  fields  in  the 
data  set.  The  increase  in  the  correlation  coefficient  after 
imposing  a 30  pixel  field  size  constraint  is  hypothesized  to 
be  due  to  the  effect  of  reducing  the  number  of  edge  or  border 
pixels  thereby  reducing  their  effect  on  the  VI  statistical 
mean  for  each  field.  Though  analysts  purposely  tried  to  avoid 
defining  fields  containing  border  pixels,  in  same  cases  the 
small  size  of  some  fields  prevented  the  definitions  of  fields 
containing  entirely  pure  wheat  pixels.  The  relationship  bet- 
ween yield  and  the  Vi's  was  relatively  high  and  stable  during 
heading. 

The  effect  of  increasing  the  field  size  constraint  during 
flowering  was  similar  to  that  found  at  heading,  although  the 
relationship  between  yield  and  the  Vi's  was  weak.  The  average 
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correlation  coefficient  was  .21,  .39,  and  .33  for  the  30,  40 
and  50  pixel  field  size  constraints,  respectively.  This 
compares  with  the  average  correlation  coefficient  at  flowering 
of  .17  based  on  all  the  fields  in  the  data  set.  Similiarly, 
as  found  during  heading,  reducing  the  effect  of  border  pixels 
improved  the  relationship  between  yield  and  the  Vi's. 

2. 2. 2. 4 Effect  of  APU  Location.  The  effect  of  segment  location  on  the 
relationship  between  yield  and  the  Vi's  was  also  investigated. 
The  field  level  data  was  stratified  according  to  its  location 
with  respect  to  the  agrophysical  units  (APU)  identified  in 
Montana  and  North  Dakota  (Figure  1-1).  The  purpose  of  comput- 
ing the  correlation  coefficients  by  APU  were  to  reduce  the 
effects  of  soils  and  climate  on  the  relationship  between  the 
Vi's  and  yield.  There  were  five  APU's  identified  in  the  two 
states  including  APU  19,  20,  21,  22  and  104.  APU  104  was  not 
included  in  the  analysis  due  to  the  low  segment  acquisition 
rate  experienced  wi thin  this  region  in  Montana. 


It  was  originally  intended  to  compute  correlation  tables, 
independently  for  each  APU  at  the  stem  extension,  heading 
and  flowering  growth  stage  intervals.  Again,  the  low  per- 
centage of  data  at  these  growth  stages  prevented  a complete 
and  thorough  analysis.  Correlations  were  only  computed  at 
flowering  for  APU  19  and  20,  at  heading  and  flowering  for  APU 
21,  and  at  stem  extension  and  flowering  for  APU  22  (refer  to 
Appendix  D for  Correlation  Tables).  A comparison  of  the  corre- 
lation coefficients  at  the  APU  level  to  those  computed  without 
respect  to  stratification,  provided  the  basis  for  analysis. 
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Only  one  APU  contained  sufficient  data  to  support  correlation 
analysis  at  stem  extension.  The  average  correlation  coefficient 
at  stem  extension  in  APU  22  was  .67,  compared  to  .36  for  the 
entire  data  set,  an  increase  of  86  percent.  Definite  conclu- 
sions cannot  be  identified  due  to  the  lack  of  sufficient  data 
at  each  of  the  other  APU's. 


Only  two  APU's  contained  a sufficient  number  of  observations 
at  heading  to  conduct  the 
coefficients  at  heading  i: 

.52,  respectively.  This  compares  with  an  average  correlation 
coefficient  at  heading  of  .58  for  the  entire  data  set.  In  this 
case,  APU  location  had  a slight  negative  affect  on  the  relation- 
ship between  yield  and  the  Vi's.  Collectively,  the  observations 
at  heading  in  the  other  APU's  supported  a strong  correlation 
between  yield  and  the  Vi's. 

Three  APU's  contained  a sufficient  number  of  observations 


analysis.  .The  average  correlation 
APU's  21  and  22  were  .45  and 


3 / 

APU  22  contained  17  observations  during  stem  extension 
4/ 

APU  21  contained  21  observations,  while  APU  22  contained  23 
observations . 
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at  flowering  to  conduct  the  analysis.  The  avenge  cgj/— 
re^tion  coefficients  at  flowering  in  APU's  19  , 20  and 

21  were  .45,  .13  and  -.18,  respectively.  This  compares 
with  an  average  correlation  coefficient  at  flowering  of 
.17  for  the  entire  data  set.  Although  the  relationship 
between  yield  and  the  Vi's  is  relatively  weak  at  flowering, 
the  stratification  by  APU  did  little  to  strengthen  this 
relationship.  The  variation  between  yield  and  the  Vi's 
was  so  great  during  this  period  that  the  relationship  was 
negative  in  APU  21  and  positive  in  APU's  19  and  20. 

Definite  conclusions  cannot  be  drawn  by  this  analysis  due 
to  the  low  volume  of  data  collected  in  each  of  the  APU's 
during  the  stem  extension,  heading  and  flowering  growth 
stage  intervals.  There  is  an  indication  that  stratification 
of  the  data  by  APU  has  little  effect  on  the  relationship 
between  yield  and  the  Vi's.  The  relationship  remained 
strongest  at  heading  in  each  APU. 

2.2.3  Yield  and  Soil  Moisture.  Surface  and  subsurface  soil  rroisture 
estimates  were  computed  frcm  the  Two  Layer  Soil  Moisture  Model. 
This  section  presents  the  correlation  analysis  results  between 
surface  and  subsurface  soil  moisture,  the  Vi's,  and  the  ground 
observed  data. 


Subsurface  soil  moisture  shewed  a stronger  relationship  to  yield 
than  was  found  with  surface  soil  moisture.  Figures  2-9  through 
2-15  shew  the  relationship  between  subsurface  soil  moisture  and 
yield.  Hie  relationship  between  yield  and  subsurface  soil  mois- 
ture peaked  at  heading  producing  a correlation  coefficient  of  .59. 
The  relationship  at  other  growth  stage  intervals  was  relatively 
strong  and  averaged  above  .45  except  at  stem  extension  when  it 
dropped  to  .15.  (Refer  to  Tables  2-1  through  2-7.)  The  soil 
moisture  estimates  frcm  the  Two  Layer  Soil  Moisture  Model  corre- 
lated reasonably  well  to  yield. 

2.2.4  Plant  Density,  Plant  Height  and  Vi's.  Tillering  and  stem  extension 
were  the  only  growth  stage  intervals  in  which  meaningful  relation- 
ships between  plant  height,  plant  density  and  each  of  the  Vi's  were 
possible.  Therefore,  the  analysis  was  limited  to  these  growth  stage 
intervals.  The  methods  of  recording  the  ground  observed  data  was 
the  principle  obstacle  to  this  analysis.  After  stem  extension  plant 
density  reached  its  peak  of  code  5 indicating  a range  of  80  to  100 
percent  ground  cover,  and  did  not  vary  during  the  remainder  of  the 
growing  season.  A similar  situation  was  found  wi th  plant  height, 
reaching  its  peak  after  stem  extension  and  not  varying  during  the 
remainder  of  the  growing  season. 


APU  19,  20  and  21  contained  18,  21  and  16  observations,  respectively. 
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FIGURE  2-Q  YIELD  VS.  SUBSURFACE  SOIL  MOISTURE  AT  PLANTING 
PLOT  OF  Y I ELD®SO I LMQ?  LEGEND:  A = 1 OBS»  B = 2 OBS*  ETC. 
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FIGURE  2-10  YIELD  VS.  SUBSURFACE  SOIL  MOISTURE  AT  TILLERING 

PLOT  OF  YIEL0*S0ILM02  LEGEND:  A = 1 OBS.  B = 2 OBS.  ETC. 
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FIGURE  2-11  YIELD  VS.  SUBSURFACE  SOI.L  MOISTURE  AT  STEM  EXTENSION 

PLOT  OF  Y I ELD*SO I LM02  LEGEND:  A = 1 ORS.  B = 2 OBS,  ETC. 
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FIGURE  2-12  YIELD  VS  SUBSURFACE  SOIL  MOISTURE  AT  HEADING 
PLOT  OF  Y I EL  D*SO I LM02  LEGEND:  A = 1 08S.  B = 2 OBS.  ETC. 
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pFJGURE  2-15  YIELD  VS.  SUBSURFACE  SOIL  MOISTURE  Al  HARVEST 
PLOT  or  Y I ELD*SO I LM02  LEGEND:  A = 1 OBS,  B = 2 OBS,  ETC. 
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2.2.4. 1 Plant  Density  and  Vi's.  The  relationship  between  plant 

density  and  the  Vi's  was  strong  at  both  tillering  and  stem 
extension.  The  average  correlation  coefficients  between 
the  Vi's  and  plant  density  were  .49  and  .67  at  tillering 
and  stem  extension,  respectively  (Tables  2-2  and  2-3). 

The  correlation  coefficients  ranged  frcm  a high  of  .51  for 
the  KVI  to  a low  of  .45  for  the  LAI  at  tillering  and  from  a 
high  of  .71  for  the  TVT6  to  a low  of  .65  for  the  DVI  at  stem 
extension.  The  VI ' s measured  plant  density  similarly. 


2. 2. 4. 2 Plant  Height  and  Vi's.  The  relationship  between  plant  height 
and  the  Vi's  was  also  strong  at  tillering  and  stem  extension. 

The  average  correlation  coefficient  between  the  Vi's  and  plant 
height  were  .66  and  .81  at  tillering  and  stan  extension,  respec- 
tively (Tables  2-2  and  2-3).  The  correlation  coefficients  at 
tillering  ranged  frcm  a high  of  .69  for  the  KVT  to  a low  of  .63 
for  the  TVI7,  and  at  stem  extension  ranged  frcm  a high  of  .84 
for  the  DVI  to  a low  of  .76  for  the  LAI.  The  Vi's  measured 
plant  height  similarly. 
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PART  3.0  REGRESSION  ANALYSIS  RESULTS 


3.1  OVERVIEW 

The  purpose  of  the  multiple  regression  analysis  was  to  identify 
functional  relationships  between  a set  of  independent  variables 
including  the  Vi's,  soil  moisture  and  a number  of  interactive  terms 
and  the  dependent  variables  including  plant  density,  plant  height  and 
yield  . This  analysis  identified  the  times  in  the  growing  year 
when  significant  relationships  to  plant  density,  plant  height  and 
yield  occurred.  Additionally,  the  analysis  identified  those  variables 
that  are  functionally  related  to  the  set  of  dependent  variables  and 
the  relative  strength  of  these  relationships. 

A stepwise  regression  procedure  was  implemented  for  each  VI  at  each 
of  the  seven  growth  stage  intervals.  A table  listing  the  independent 
variables  used  in  the  regression  analysis  is  presented  as  Appendix  E. 

Hie  Vi's  were  never  jointly  analyzed  in  a regression,  but  rather  were 
individually  regressed  against  the  set  of  dependent  variables.  The 
results  of  this  analysis  are  presented  in  the  following  sections. 

3.2  RELATIONSHIPS  TO  YIELD 

3.2.1  Plant ing.  The  VI ' s did  not  enter  into  the  regression  relationship 
defined  at  planting  (Table  3-1).  Instead,  subsurface  soil  moisture 
raised  to  the  fourth  power  was  the  only  significant  variable 
entering  into  the  regression  relationship.  A coefficient  of 
determination  of  .32  and  a standard  deviation  of  8.12  indicated 

a relatively  weak  relationship  to  yield.  The  relative  strength 
of  the  relationship  at  planting  was  expected,  because  it  is  too 
early  in  the  growing  season  to  accurately  explain  final  yield. 
Although  the  relationship  was  weak  at  planting,  it  was  expected 
that  subsurface  soil  moisture  would  have  some  relationship  to 
final  yield  at  planting. 

3.2.2  Ti 1 ler ing.  The  strength  of  the  relationship  at  tillering  was  also 
weak  and  produced  an  average  coefficient  of  determination  and 
standard  deviation  for  all  the  Vi's  of  .28  and  8.40,  respectively 
(Table  3-2).  The  independent  variables  entering  consistently 
into  the  regression  equations  were  a dimmy  variable  differ- 
entiating between  spring  and  winter  wheat,  subsurface  soil 
moisture  raised  to  the  third  power,  the  VI  multiplied  by  the  sun  of 
surface  and  subsurface  soil  moisture,  and  the  VI  added  to  subsurface 


7 J 

Since  there  was  a perfect  correlation  between  the  KVI  and  the  GVI  (due 
to  sun  angle  and  haze  correction)  as  reported  in  Part  2.0,  the  GVI 
was  deleted  from  the  list  of  Vi's  analyzed  by  regression  analysis. 

The  results  for  both  Vi's  would  have  been  similar. 

g J 

A dummy  variable  is  used  to  measure  the  difference  between  factors 
that  cannot  easily  be  quantified.  In  this  case  the  factor  being 
measured  is  the  effect  of  wheat  type  on  final  yield. 
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TABLE  3-1  YIELD  REGRESSION  RELATIONSHIP  AND  ASSOCIATED  STATISTICS 
FOR  EACH  VI  AT  PLANTING 


VI 

1 

N 

2 

R2 

3 

STD  DEV 

4 

B-VALUE 

5 

STD  ERROR 

6 

F-VALUE 

7 

PROB  > F 

ALL  Vi's 

88 

.32 

8.  12 

- 

- 

40.22 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

21.862 

- 

— 

- 

SM2F0 

0.002 

.0004 

40.22 

.0001 

lMNM  is  the  number  of  observations. 

2"R2"  is  the  coefficient  of  determination  indicating  the  amount  of  variation 
in  the  dependent  variable  (yield)  explained  by  the  independent  variable(s). 

"R2"  ranges  from  0.00  to  1.00, 

""STD  DEV"  is  the  average  deviation  between  the  observed  and  estimated  yield. 

^"B-VALUES"  are  the  coefficients  or  weights  applied  to  each  of  the  independent 
variables.  Additionally,  the  intercept  is  also  listed  under  "B-VALUE." 

5"STD  ERROR"  is  the  standard  error  of  the  B-Values. 

6"F-VALUE"  is  used  to  test  whether  the  linear  relationship  is  significant.  Addi- 
tionally, it  is  used  to  test  whether  a variable  meets  the  significance  criteria 
and  should  be  included  in  the  model  (partial  F-Test) . The  F-Value  is  the  ratio 
of  the  regression  mean  square  to  the  error  mean  square. 

7"PR0B  >F"  is  the  significance  probability  of  the  F-Value. 
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TABLE  3-2  YIELD  REGRESSION  RELATIONSHIP  AND  ASSOCIATED  STATISTICS 
FOR  EACH  VI  AT  TILLERING 


VI 

N 

r2~ 

STD  DEV 

B-VALUE 

STD  ERROR 

F-VA1UE 

PROB  >F 

AVI 

196 

. 24 

8.66 

_ 

- 

19.97 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

33.46 

- 

- 

- 

Whea  t 

- 

- 

- 

-3.80 

1.452 

6.86 

.0095 

SMICU 

- 

- 

- 

-6.54 

2.258 

8.47 

.0040 

SM2CU 

- 

- 

- 

0.02 

0.003 

4 5.06 

.0001 

DVI 

196 

.29 

8.35 

- 

- 

20.06 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

19.61 

- 

- 

- 

Wheat 

- 

- 

- 

-3.79 

1.430 

7.04 

.0086 

SM2CU 

- 

- 

- 

0.02 

0.003 

49.61 

.0001 

DVIMSM3 

- 

- 

- 

-0.25 

0.056 

21.43 

.0001 

DVIASM2 

- 

- 

- 

1.82 

0.426 

18.  20 

.0001 

LAI 

196 

.30 

8.30 

- 

- 

20.89 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

18.00 

- 

- 

- 

Wheat 

- 

- 

- 

-3.80 

1.450 

6.87 

.0095 

SM2CU 

- 

- 

- 

0.02 

0.003 

27.53 

.0001 

LAIMSM3 

- 

- 

- 

-0.32 

0.  066 

23.52 

.0001 

LATASM2 

- 

- 

- 

2.26 

0.  511 

19.64 

.0001 

KVI 

196 

.30 

8.33 

- 

- 

20.26 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

14.33 

- 

- 

- 

Wheat 

- 

- 

- 

-3.80 

1.429 

7.08 

.0084 

SM2CU 

- 

- 

- 

0.02 

0.003 

36.67 

.0001 

K\'  IMSM3 

- 

- 

- 

-0.39 

0.082 

22.25 

.0001 

KV  I ASM 2 

- 

- 

- 

2.73 

0.628 

18.95 

.0001 

TVI6 

196 

.27 

8.49 

- 

- 

23.44 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

19.35 

- 

- 

- 

Wheat 

- 

- 

- • T - 

-2.97 

1.38 

4.62 

.0328 

S0ILM02 

- 

- 

- 

8.60 

1.27 

46.20 

.0001 

TVIMSM3 

- 

- 

- 

-5  25 

1.24 

18.03 

.0001 

TV  17 

196 

.25 

' 8.57 

__  i 

21.69 

.0001 

VARIABLES 

In  tercept 

- 

- 

- 

10.79 

- 

- 

- 

Wheat 

- 

- 

- 

-2.37 

1.40 

2.87 

.0919 

TVLKSM3 

- 

- 

- 

-4.81 

1.23 

15.26 

.0001 

TV1ASM2 

- 

- 

- 

7.21 

1.04 

47.95 

.0001 

PVI6 

196 

.31 

8.24 

- 

- 

21.81 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

4.00 

- 

- 

- 

Wheat 

- 

- 

- 

-4.12 

1.38 

8.98 

.0031 

PV  LMSM3 

- 

- 

- 

-0.39 

0.  09 

17.09 

.0001 

PVIASM1 

- 

- 

- 

-2.86 

0.53 

28.99 

.0001 

FVIASM2 

- 

- 

- 

5.53 

0.69 

63.46 

.0001 

PVI7 

196 

.30 

8.28 

- 

- 

27.76 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

0.39 

- 

- 

- 

Wheat 

- 

- 

- 

-3.61 

1.39 

6.71 

.0103 

PVLMSM3 

- 

- 

- 

-0.81 

0.09 

79.49 

.0001 

* PVIASM2 

- 

- 

- 

5.68 

0.  70 

65.42 

.0001 
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soil  iroisture.  There  was  minimal  variation  in  the  variables 
entering  into  the  yield  relationship  when  using  each  VI.  The 
similarity  of  the  regression  relationships  for  each  VI  was  ex- 
pected due  to  the  high  correlation  among  the  Vi's.  Again,  soil 
moisture  proved  to  be  an  important  variable  in  explaining  the 

variation  in  yield.  Hie  interaction  between  the  Vi's  and  soil 
moisture  also  proved  to  be  extrenely  significant  in  explaining 
the  variation  in  wheat  yield. 

Historically,  winter  wheat  because  of  its  longer  growing  season 
along  with  physiologically  based  differences  produces  a higher 
yield  when  compared  to  the  yields  produced  by  spring  wheat.  The 
regression  relationship  defined  at  tillering  depicted  this  relation- 
ship. A code  variable  was  defined  to  differentiate  between  wheat 
types  with  a value  of  one  given  to  the  winter  wheat  fields  and  a 
value  of  two  given  to  the  spring  wheat  fields.  The  negative 
coefficient  preceding  the  wheat  type  variable  indicates  that  at  a 
specific  growth  stage  given  a value  of  a VI,  the  yield  of  a field 
planted  to  spring  wheat  will  be  less  than  the  yield  planted  to 
a field  of  winter  wheat. 

3.2.3  Stan  Extension.  The  relationship  between  yield  and  the  independent 
variables  improved  at  stem  extension.  The  coefficient  of  deter- 
mination and  the  standard  deviation  for  all  the  Vi's  averaged  .41 
and  5.87,  respectively  (Table  3-3).  The  standard  deviation  ranged 
from  a low  of  5.51  for  the  LAI  to  a high  of  6.32  for  the  TVI6. 

The  independent  variables  consistently  entering  into  the  regression 
equation  for  each  VI  were  a dirrmy  variable  indicating  the  wheat 
type  (i.e.,  winter  or  spring  wheat),  the  VI  multiplied  by  the  sun 
of  surface  and  subsurface  soil  moisture,  and  the  VI  multiplied 
by  the  growth  stage  within  the  sten  extension  interval.  The  LAI 
multiplied  by  subsurface  soil  moisture  was  the  most  significant 
variable  found  in  the  relationships  defined  for  the  TVI6  and  TVI7. 

The  VI  multiplied  by  the  growth  stage  was  not  significant  in  the 
case  of  the  yield  regression  relationship  defined  for  the  LAI. 

During  stem  extension  the  wheat  plant  is  undergoing  a significant 
amount  of  growth  (Figure  1-2).  Plant  density  within  the  fields 
is  increasing  at  a maximum  rate.  The  VI  multiplied  by  the  growth 
stage  is  understandably  an  important  variable  in  explaining  the 
variation  in  yield.  This  term  adjusts  the  VI  for  the  growth  stage 
within  the  stem  extension  interval.  For  example,  a VI  of  25  at 
growth  stage  6.0  would  indicate  a greater  potential  yield  when 
compared  to  the  same  VI  at  growth  stage  8.0. 

The  interaction  between  the  VI  and  surface  and  subsurface  soil 
moisture  also  proved  to  be  significant.  In  fact,  in  most  of  the 
relationships  this  term  was  the  must  significant  variable  to  yield 
(refer  to  significance  levels  for  each  of  the  independent  variables). 

The  dirrmy  variable  differentiating  between  wheat  types  proved  to 
be  significant  at  step  extension,  as  was  the  case  at  tillering. 
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TABLE  3-3  YIELD  REGRESSION  RELATIONSHIP  AND  ASSOCIATED  STATISTICS 
FOR  EACH  VI  AT  STEM  EXTENSION 


VI 

N 

R2 

STD  DEV 

B- VALUE 

STD  ERROR 

F-VAIUE 

PROB  >F 

AVI 

31 

.46 

3.68 

_ 

_ 

7.68 

.0007 

VARIABLES 

Intercept 

- 

- 

- 

39.57 

- 

- 

- 

Wheat 

- 

- 

- 

-8.32 

2.53 

10.84 

.0028 

AVIMSM3 

- 

- 

- 

0. 14 

0.04 

12.53 

.0015 

AV1MGRO 

- 

- 

- 

-0.08 

0.03 

7.69 

.0100 

DVI 

31 

.44 

5.79 

- 

- 

7.10 

.0011 

VARIABLES 

Intercept 

- 

- 

- 

36.78 

- 

- 

- 

Wheat 

- 

- 

- 

-8.33 

2.59 

10.35 

.0034 

DV1MSM3 

- 

- 

- 

0. 11 

0.03 

11.64 

.0021 

DVIMGRO 

- 

- 

- 

-0.  53 

0.02 

6.40 

.0176 

LAI 

31 

.49 

5.51 

- 

- 

8.78 

.0003 

VARIABLES 

Intercept 

- 

- 

- 

35.71 

- 

- 

- 

Wheat 

- 

- 

- 

-8.98 

2.53 

12.59 

.0015 

LAIMSM2 

- 

- 

- 

0.17 

0.04 

15.82 

.0005 

LAIMGRO 

- 

- 

- 

-0.07 

0.02 

8.40 

.0074 

KVI 

31 

.43 

5.82 

- 

- 

6.90 

.0013 

VARIABLES 

Intercept 

- 

- 

- 

35.12 

- 

- 

- 

Wheat 

- 

- 

- 

-8.34 

2.63 

10.08 

. 0039 

KVLMSM3 

- 

- 

- 

0.16 

0.05 

11.56 

.0021 

KVI MG RO 

- 

- 

- 

-0.  07 

0.03 

6.20 

.0193 

TVI6 

31 

.31 

6.32 

- 

- 

6.26 

.0057 

VARIABLES 

Intercept 

- 

- 

- 

10.46 

- 

- 

- 

Wheat 

- 

- 

- 

-6.33 

2.60 

5.94 

.0214 

TVI6MSM3 

- 

- 

- 

3.88 

1.44 

7.23 

.0119 

TV  17 

31 

.37  * 

* 6.02- 

- 

' - 

• 8.33 

.0014 

VARIABLES 

Intercept 

- 

- 

- 

12.02 

- 

- 

- 

Wheat 

- 

- 

- 

• -5.98 

2.47 

5.86 

.0222 

TVI7SM3 

- 

- 

- 

4.31 

1.31 

10.84 

.0027 

PVI6 

31 

.38 

6.07 

- 

- 

5.64 

.0039 

VARIABLES 

Intercept 

- 

- 

- 

35.06 

- 

- 

- 

Wheat 

- 

- 

- 

-8.36 

2.76 

9.18 

.0053 

PVI6MSM3 

- 

- 

- 

0.17 

0.06 

9.20 

.0053 

PVI6MGRO 

- 

- 

- 

-0.08 

0.04 

4.94 

.0348 

PVI7 

31 

.44 

5.79 

- 

- 

7.09 

.0012 

VARIABLES 

Intercept 

- 

- 

- 

36.70 

- 

- 

- 

Wheat 

- 

- 

- 

-8.31 

2.59 

10.31 

.0034 

PVI7MSM3 

- 

- 

- 

0.28 

0.08 

11.63 

.0021 

PVI7MGRO 

- 

- 

- 

-0. 14 

0.05 

6.  39 

.0176 

> 
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3.2.4  Heading.  The  relationship  between  yield  and  the  independent  vari- 
ables peaked  at  heading.  The  coefficient  of  determination  and  the 
standard  deviation  between  the  observed  and  estimated  yields 
averaged  .57  and  5.75,  respectively  (Table  3-4).  The  standard 
deviation  ranged  from  a low  of  5.50  for  the  LAI  to  a high  of  5.90 
for  the  PVI6.  For  each  VI,  there  was  very  little  variation  in  the 
variables  entering  the  regression  relationship  and  the  relative 
strength  of  the  relationship  at  heading.  This  was  due  to  the 
relatively  high  correlation  among  the  Vi's  at  heading. 

The  independent  variable  consistently  entering  into  the  regression 
equation  for  each  VI  and  the  variable  which  proved  to  be  the  most 
significant  was  the  VI  multiplied  by  subsurface  soil  moisture.  The 
durrmy  variable  differentiating  between  winter  and  spring  wheat  was 
also  significant  in  the  AVI,  DVI  and  PVI7  yield  relationships.  The 
VI  added  to  subsurface  soil  moisture  also  proved  to  be  significant 
in  the  TVI6  and  TVI7  relationships. 

Figure  3-1  shows  a ccmparison  between  the  observed  yields  and  the 
yields  estimated  by  the  LAI  derived  yield  relationship.  The  term 
"LAIMSM2"  (LAI  multiplied  by  subsurface  soil  moisture)  proved  to  be 
the  only  significant  variable  in  the  yield  relationship  defined  at 
heading.  The  difference  between  the  observed  and  estimated  yields 
averaged  5.50  bushels /acre. 

A significant  improvement  in  the  correlation  coefficient  at  heading 
between  the  LAI  and  yield  resulted  frcm  an  increasing  field  size 
limitation  (Section  2. 2. 2. 3).  This  prcrrpted  further  analysis 
investigating  the  effect  of  field  size  on  the  yield  regression 
relat ionship. 

The  minimum  field  size  limitations  investigated  were  established 
at  30,  35,  40  and  45  pixels.  For  example,  a 30  pixel  minimum 
field  size  limitation  prevented  any  fields  less  than  30  pixels  in 
size  frcm  entering  the  yield  regression  analysis. 

The  analyses  identified  the  "LAIMSM3"  (LAI  multiplied  by  the  sum 
of  surface  and  subsurface  soil  moisture)  as  the  only  significant 
variable  explaining  the  variation  in  yield  as  a result  of  imposing 
the  minimum  field  size  limitations.  The  regression  relationship 
improved  significantly  by  implementing  the  field  size  constraints. 
The  coefficient  of  determinat ions  were  .81,  .80,  .80  and  .78  by 
imposing  the  30,  35,  40  and  45  pixel  field  size  constraints, 
respectively,  compared  to  a .61  for  all  the  field  data.  Hie 
standard  deviations  between  the  observed  and  estimated  yields  were 
4.08,  4.18,  4.12  and  4.37  the  30,  35,  40  and  45  pixel  field  size 
constraints,  respectively,  compared  to  a 4.80  for  all  the  field 
data.  Figure  3-2  shows  a comparison  between  the  observed  yields 
and  the  yields  estimated  by  the  LM  derived  relationship  for 
fields  greater  than  30  pixels  insize. 

The  number  of  observations  used  to  build  the  regression  relation- 
ship after  imposing  the  30,  35  40  and  45  pixel  field  size  con- 
straints were  29,  26,  21  and  18  respectively,  compared  to  55 
observations  for  the  entire  data  set  at  heading.  Over  45  percent 
of  the  fields  in  the  data  set  were  less  than  30  pixels  in  size. 
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TABLE  3'4  YIELD  REGRESSION  RELATIONSHIP  AND  ASSOCIATED  STATISTICS 
FOR  EACH  VI  AT  HEADING 


VI 

N 

R2 

STD  DEV 

B-VALUE 

STD  ERROR 

F -VALUE 

PROB  >F 

AVI 

53 

.57 

5.81 

- 

_ 

34.41 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

25.33 

- 

- 

- 

Wheat 

- 

- 

- 

-4.19 

1.72 

5.94 

.0183 

AVIMSM2 

— 

0.06 

0.01 

45.43 

.0001 

DVI 

55 

.56 

5.84 

_ 

_ 

33.71 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

23.55 

- 

- 

- 

Wheat 

- 

- 

- 

-3.86 

1.75 

4.90 

.0313 

DVIMSM2 

“ 

- 

0.05 

0.01 

44.31 

.0001 

LAI 

55 

.61 

5.50 

_ 

81.68 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

15.86 

- 

- 

- 

LAI MSM2 

"" 

0.08 

0.01 

81.68 

.0001 

KVI 

55 

.55 

5.83 

- 

- 

66.06 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

15.27 

- 

- 

- 

KVIMSM2 

— 

0.09 

0.01 

66.06 

.0001 

TVI6 

55 

.60 

5.59 

- 

- 

39.19 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

36.94 

- 

- 

- 

TVIMSM2 

- 

- 

- 

17.61 

3.13 

31.59 

.0001 

TVIASM2 

— 

-18. 52 

3.79 

23.93 

.0001 

TVI7 

55 

.58 

5.73 

- 

- 

36.11 

.0001 

VARIABLES 

Intercept 

23.16 

TVIMSM2 

- 

- 

- 

10.18 

1.97 

26.76 

.0001 

TVIASM2 

- 

- 

-8.03 

2.10 

14.57 

.0004 

PVI6 

55 

.55 

5.90 

- 

- 

64.03 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

15.09 

- 

- 

- 

AVIMSM2 

“ 

0.10 

.01 

O' 

O 

U> 

.0001 

PVI7 

55 

.56 

5.84 

- 

- 

33.71 

.0001 

VARIABLES 

Intercept 

23.56 

Wheat 

- 

- 

- 

-3.86 

1.75 

4.90 

.0313 

AVIMSM2 

- 

- 

- 

0.14 

0.02 

44.31 

.0001 

« 
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FIGURE  3-2  OBSERVED  AND  ESTIMATED  YIELDS 
vs.  LAIMSM3  AT  HEADING,  FIELDS  CREATER  THAN  30  riXELS 
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The  improvement  in  the  yield  regression  relationship  at  heading 
is  believed  to  be  due  to  the  reduction  in  the  effect  border  or 
edge  pixels  have  on  the  VI  statistical  mean  at  the  field  level. 

That  is,  as  the  minirrum  field  size  is  increased,  the  percentage 
of  edge  pixels  to  pure  pixels  is  reduced.  The  analysis  indicated 

that  the  yield  relationship  inproved  significantly  after  imposing 
the  30  pixel  field  size  constraint,  and  remained  stable  as  the 
field  size  constraint  was  increased  to  35,  40  and  45  pixels. 

Further  constraints  could  not  be  imposed  due  to  an  insufficient 
number  of  field  observations  greater  than  45  pixels  in  size. 

3.2.5  Flowering.  The  relationship  between  yield  and  the  independent 
variables  was  also  relatively  strong  at  flowering.  The  coefficient 
of  determination  and  the  standard  deviation  between  the  observed 

and  estimated  yields  averaged  .62  and  6.00,  respectively  (Table  3-5). 
The  standard  deviation  ranged  frcm  a low  of  5.91  for  the  TVI6  to  a 
high  of  6.06  for  the  PVI6.  There  was  a relatively  small  difference 
between  the  strength  of  the  yield  relationship  defined  at  heading 
and  flowering.  In  fact,  the  difference  in  the  average  standard 
deviation  at  heading  and  flowering  was  0.25,  a difference  of  4 per- 
cent. Again,  there  was  little  variation  in  the  independent  variables 
entering  the  regression  relationship  and  the  relative  strength  of 
the  relationship  for  each  VI  due  to  the  high  correlation  among  the 
Vi's  at  flowering. 

The  sum  of  surface  and  subsurface  soil  moisture  raised  to  the 
fourth  power  proved  to  be  the  most  significant  variable  entering 
into  the  yield  relationship  for  each  VI.  Other  variables  which 
proved  to  be  significant  were  the  VI  multiplied  by  surface  soil 
moisture,  the  dummy  variable  differentiating  between  wheat  types 
and  other  soil  moisture  derived  variables. 

3.2.6  Ripening.  At  ripening,  the  relationship  between  yield  and  the 
independent  variables  was  weakest  and  produced  an  average  co- 
efficient of  determination  and  standard  deviation  between  the 
observed  and  estimated  yields  of  .26  and  9.52,  respectively 
(Table  3-6).  The  standard  deviation  ranged  frcm  a low  of  9.45 
for  the  KVI  to  a high  of  9.70  for  the  PVI7.  There  was  little 
variation  in  the  independent  variables  entering  the  regression 
relationship  and  the  relative  strength  of  the  relationship  for 
each  VI  at  ripening. 

Variables  defined  frcm  subsurface  soil  moisture  proved  to  be  the 
most  significant  to  yield.  Those  variables  included  the  sun  of 
surface  and  subsurface  soil  moisture  raised  to  the  fourth  power, 
and  the  VI  multiplied  by  subsurface  soil  moisture  and  subsurface 
soil  moisture.  Other  variables  which  were  significant  included 
the  code  variable  differentiating  between  wheat  types  and 
the  VI  multiplied  by  growth  stage. 

3.2.7  Harvest . The  relationship  between  yield  and  the  independent 
variables  was  stronger  at  harvest  then  was  found  at  ripening. 

The  coefficient  of  determination  and  the  standard  deviation 
between  the  observed  and  estimated  yields  averaged  .42  and  7.36, 
respectively  (Table  3-7).  The  standard  deviation  ranged  frcm 
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TABLK  3~5  YIELD  REGRESSION  RFLATIONSHI  P AND  ASSOCIATED  STATISTICS 
FOR  EACH  VI  AT  FLOWERING 


VI 

N 

R2 

STD  DEV 

B-VALUE 

STD  ERROR 

F-VALUE 

PKOB  >F 

AVI 

99 

.61 

6.04 

- 

- 

29.69 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

26.00 

- 

- 

- 

SM2SQ 

- 

- 

- 

-4.18 

1.02 

16.78 

.0001 

SM2FO 

- 

- 

- 

0.07 

0.01 

27.18 

.0001 

SM3CU 

- 

- 

- 

1.43 

0.27 

28.20 

.0001 

SM3FO 

- 

- 

-0.19 

0.03 

34.51 

.0001 

AVIMSMl 

- 

- 

- 

0.  48 

0.25 

3.49 

.0648 

DVI 

99 

.61 

6.04 

- 

- 

29.72 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

26.00 

- 

- 

- 

SM2SQ 

- 

- 

- 

-4.18 

1.02 

16.85 

.0001 

SM2FO 

- 

- 

- 

0.07 

0.01 

26.74 

.0001 

SM3CU 

- 

- 

- 

1.43 

0.27 

28.29 

.0001 

SM3FO 

- 

- 

- 

-0.18 

0.03 

34.66 

.0001 

DVIMSM1 

- 

- 

- 

0.40 

0.21 

3.55 

.0627 

LAI 

99 

. 64 

5.97 

- 

- 

27.31 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

62.02 

- 

- 

- 

Wheat 

- 

- 

- 

-5.31 

2.57 

4.27 

.0415 

S0ILM02 

- 

- 

- 

-29.60 

5.45 

29.51 

.0001 

SM2SQ 

- 

- 

- 

5.90 

0.93 

40.52 

.0001 

SM3FO 

- 

- 

- 

-0.03 

0.01 

50.73 

.0001 

LATMSMl 

- 

- 

- 

0.91 

0.19 

22.80 

.0001 

LAIMSM3 

- 

- 

- 

0.04 

0.01 

6.57 

.0120 

KVI 

99 

.61 

6.05 

- 

- 

29.50 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

26.00 

- 

- 

- 

SM2SQ 

- 

- 

- 

-4.19 

1.02 

16.84 

.0001 

SM2FO 

- 

- 

- 

0.08 

0.01 

24.96 

.0001 

SM3CU 

- 

- 

- 

1.43 

0.27 

28.27 

.0001 

SM3F0 

- 

- 

- 

-0. 19 

0.03 

34.09 

.0001 

KVLMSM1 

- 

- 

- 

0.66-  J 

- 0.38 

3.10 

.0816 

TVI6 

99 

.63 

5.91 

- 

- 

26.68 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

114.04 

- 

- 

- 

Wheat 

- 

- 

- 

-4.82 

2.56 

3.54 

.0630 

SM1FO 

- 

- 

- 

32.30 

6.44 

25.17 

.0001 

SM2SQ 

- 

- 

- 

5.79 

0.91 

39.98 

.0001 

SM3FO 

- 

- 

- 

-0.04 

0.01 

50.38 

.0001 

TVLMSM2 

- 

- 

- 

5.85 

3.65 

18.80 

.0001 

TVIASM2 

- 

- 

- 

-46.40 

7.63 

36.96 

.0001 

TVI7 

99 

.63 

5.95 

- 

- 

26.19 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

102.45 

- 

- 

- 

Wheat 

- 

- 

- 

-4.83 

2.58 

3.52 

.0639 

SM2SQ 

- 

- 

- 

6.16 

0.96 

40.91 

.0001 

SM3FO 

- 

- 

- 

-0.04 

0.01 

50.60 

.0001 

TVIMSM1 

- 

- 

- 

32.33 

6.32 

26.12 

.0001 

TVIMSM2 

- 

- 

- 

10.89 

2.36 

21.26 

.0001 

TVIASM2 

- 

- 

- 

-40.74 

6.93 

34.56 

.0001 

PVI6 

99 

.61 

6.06 

- 

- 

29.40 

.0001 

VARIABLES 

. Intercept 

- 

- 

- 

26.00 

- 

- 

- 

SM2SQ 

- 

- 

- 

-4.19 

1.02 

16.77 

.0001 

SM2FO 

- 

- 

- 

0.08 

0.01 

25.22 

.0001 

SM3CU 

- 

- 

- 

1.43 

0.27 

28.15 

.0001 

SM3F0 

- 

- 

- 

-0. 19 

0.03 

33.86 

.0001 

PVIMSM1 

- 

- 

- 

0.70 

0.41 

2.91 

.0912 

» 
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YIELD  REGRESSION  RELATIONSHIP  AND  ASSOCIATED  STATISTICS 
FOR  EACH  VI  AT  FLOWERING  (Continued) 
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TABLE  3-6  YIELD  REGRESSION  RELATIONSHIP  AND  ASSOCIATED  STATISTICS 
FOR  EACH  VI  AT  RIPENING 


VI 

N 

R2 

STD  DEV 

B-  VALUE 

STD  ERROR 

F-VALUE 

PROB  >F 

AVI 

293 

.26 

9.52 

- 

- 

17.06 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

5.76 

- 

- 

- 

Wheat 

- 

- 

- 

4.40 

1.33 

10.87 

.0001 

S01LM02 

- 

- 

- 

-10.95 

5.16 

4.50 

.0348 

SMISQ 

- 

- 

- 

-10.19 

5.03 

4.10 

.0439 

SM3FO 

- 

- 

- 

-0.004 

0.001 

11.02 

.0010 

AVI  ASM 2 

“ 

17.07 

5.06 

11.38 

.0008 

DVI 

293 

.27 

9.51 

- 

- 

14.86 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

-236.08 

- 

- 

- 

CROWSTA 

- 

- 

- 

21.81 

11.01 

3.92 

.0485 

Wheat 

- 

- 

- 

4.05 

1.32 

9.37 

.0024 

S0ILM02 

- 

- 

- 

-22.21 

8.21 

7.31 

.0073 

DVLASM1 

- 

- 

- 

-0.004 

0.001 

11.34 

.0009 

DVIASM2 

- 

-9.03 

4.10 

4.86 

.0283 

LAI 

293 

.27 

9.47 

- 

- 

15.36 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

32.10 

- 

- 

- 

LAI 

- 

- 

- 

11.07 

5.04 

4.82 

.0289 

Wheat 

- 

- 

- 

4.83 

1.36 

12.55 

.0005 

S0ILM02 

- 

- 

- 

-3.68 

3.80 

0.94 

.3333 

SM2SQ 

- 

- 

- 

1.21 

0.57 

4.51 

.0345 

SM3FO 

- 

- 

- 

-0.01 

0.004 

14.58 

.0002 

LAIMSM2 

- 

- 

0.26 

0.07 

12.86 

.0004 

KVI 

293 

.28 

9.45 

- 

- 

15.67 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

-356.71 

- 

- 

- 

KVI 

- 

- 

- 

35.53 

10.13 

12.31 

.0005 

CROWSTA 

- 

- 

- 

32.67 

12.92 

6.39 

.0120 

WTieat 

- 

- 

- 

4.03 

1.30 

9.64 

.0021 

S0ILM02 

- 

- 

- 

6.27 

0.86 

52.62 

.0001 

SM1CU 

- 

- 

- 

-12.92 

5.76 

5.03 

.0257 

SM3FO 

-0.004 

0.001 

11.62 

.0007 

TVI6 

293 

.28 

9.46 

- 

- 

13.68 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

-3810.52 

- 

- 

- 

GROWSTA 

- 

- 

- 

360.69 

142.70 

6.39 

.0120 

Wheat 

- 

- 

- 

4.54 

1.37 

10.96 

.0010 

S0ILM02 

- 

- 

- 

-3630.85 

1484.98 

5.98 

.0150 

SM2SQ 

- 

- 

- 

1.26 

0.57 

4.88 

.0280 

SM3FO 

- 

- 

- 

-0.01 

0.004 

13.15 

.0003 

TVIMSM2 

- 

- 

- 

26.52 

13.31 

3.97 

.0472 

TVIASM2 

- 

- 

3600.75 

1491.04 

5.83 

.0164 

TV17 

293 

.25 

9.60 

- 

- 

16.01 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

142.54 

- 

- 

- 

Wheat 

- 

- 

- 

4.82 

1.41 

11.66 

.0007 

S0ILM02 

- 

- 

- 

-19.18 

7.36 

6.78 

.009  7 

SM2SQ 

- 

- 

- 

0.99 

0.57 

2.97 

.0862 

SM3FO 

- 

- 

- 

-0.01 

0.004 

11.06 

.0010 

TV 1MSM2 

- 

“ 

22.40 

7.11 

9.92 

.0018 

/ 
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TABLE  3-6  YIELD  REGRESSION  RELATIONSHIP  AND  ASSOCIATED  STATISTICS 
FOR  EACH  VI  AT  RIPENING  (Continued) 


VI 

N 

R2 

STD  DEV 

B-VALUE 

STD  KRR01 

F -VALUE 

TROB  > F 

PV16 

293 

.27 

9. 47 

- 

- 

21.08 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

12.88 

- 

- 

- 

V.'Heat 

- 

- 

- 

3.35 

1.28 

6.86 

.0093 

S0ILM03 

- 

- 

- 

-6.09 

3.21 

3.61 

.0583 

SM3FO 

- 

- 

- 

-0.004 

0.001 

10.93 

.0011 

PVLASM2 

12.26 

3.11 

15.58 

.0001 

rvi7 

293 

.23 

9.70 

- 

- 

21.66 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

6.92 

- 

- 

- 

Wheat 

- 

- 

- 

3.42 

1.33 

6.64 

.0105 

S0ILH02 

- 

- 

- 

6.89 

0.92 

56.60 

.0001 

SM3F0 

- 

- 

- 

-0.004 

0.001 

11.42 

.0008 

PVIASM3 

- 

- 

- 

-0.  52 

0.14 

13.10 

.0003 

■ 
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TABLE  3"7  YIELD  REGRESSION  RELATIONSHIP  AND  ASSOCIATED  STATISTICS 
EOK  EACH  VI  AT  HARVEST 


VI 

N 

R2 

STD  DEV 

B -VALUE 

STD 

AVI 

218 

.43 

7.35 

- 

VARIABLES 

Intercept 

- 

- 

- 

24.24 

Wheat 

- 

- 

- 

-6.92 

1. 

S0I1M02 

- 

- 

7.13 

0. 

SM3FO 

- 

- 

- 

-0.01 

0. 

AV1MSM2 

- 

- 

- 

-0.  21 

0. 

DVI 

218 

.43 

7.35 

- 

VARIABLES 

Intercept 

- 

- 

- 

26.86 

Wheat 

- 

- 

- 

-8.14 

1. 

S0ILM02 

- 

- 

- 

8.50 

0. 

SM3F0 

- 

- 

- 

-0.01 

0. 

DV1MSM2 

- 

- 

- 

-0.18 

0. 

LAI 

218 

.42 

7.42 

- 

VARIABLES 

Intercept 

- 

- 

- 

28.94 

Wheat 

- 

- 

- 

-5.99 

1. 

S0IU102 

- 

- 

- 

5.11 

1. 

SM3CU 

- 

- 

- 

0.29 

0. 

SM3FO 

- 

- 

- 

-0.05 

0. 

LAIMSM2 

- 

- 

- 

-0.37 

0. 

KVI 

218 

.43 

7.32 

- 

VARIABLES 

28.36 

Intercept 

- 

- 

Wheat 

- 

- 

- 

-8.46 

1. 

S0ILM02 

- 

- 

8.79 

0. 

SM3F0 

- 

- 

- 

-0.10 

0. 

KVIMSM2 

- 

- 

- 

-0.  31 

0. 

TVI6 

218 

.42 

7.36 

- 

VARIABLES 

-59.04 

Intercept 

- 

- 

Wheat 

- 

- 

- 

-7.13 

1 

SM3FO 

- 

- 

- 

-0. 10 

0 

TVIMSM2 

- 

- 

- 

-72.49 

8 

TVIASM2 

- 

- 

- 

82. 14 

8 

TVI7 

218 

.43 

7.32 

- 

VARIABLES 

25.21 

Intercept 

- 

Wheat 

- 

- 

- 

-6.79 

1 

S0ILM02 

- 

- 

35.77 

3 

SM3FO 

- 

- 

- 

-0.01 

0 

TVIMSM2 

- 

- 

- 

-35.69 

4 

PV16 

218 

.40 

7.52 

- 

VARIABLES 

30.70 

Intercept 

- 

- 

“ 

Wheat 

- 

- 

- 

-8.86 

1 

SOILM02 

- 

- 

“ 

8.39 

0 

SM3FO 

- 

- 

- 

-0.006 

0 

PVIMSM3 

- 

- 

-0.30 

0 

PVI7 

218 

.44 

7.23 

- 

VARIABLES 

28.32 

Intercept 

- 

- 

- 

Wheat 

- 

- 

- 

-8.55 

1 

S01LM02 

- 

- 

8.48 

0 

SM3F0 

- 

- 

- 

-0. 10 

0 

PVIMSM3 

- 

- 

- 

-0.52 

0 

F -VALUE 


,48 


,44 

,55 


,48 


,46 


39.62 


22.99 
60.  47 
15.71 
70.53 


39.  78 


30.32 

81.96 

12.77 

71.03 


30.64 


17.23 

10.92 

5.32 

7.70 

59.43 


40.47 


32.57 

86.81 

9.84 

73.12 


39.42 


24.15 

8.89 

74.27 

86.55 


ROB  >F 


40.51 


22.43 

103.49 

14.22 

73.25 


35.58 


32.43 

76.17 

3.76 

58.19 


42.76 


34.16 
84.63 
9.65 
80. 12 


.0001 


.0001 

.0001 

.0001 

.0001 


.0001 


.0001 

.0001 

.0004 

.0001 


.0001 


.0001 

.0011 

.0220 

.0060 

.0001 


.0001 


.0001 

.0001 

.0020 

.0001 


.0001 


.0001 

.0032 

.0001 

.0001 


.0001 


.0001 

.0001 

.0002 

.0001 


.0001 


.0001 

.0001 

.0540 

.0001 


.0001 


.0001 

.0001 

.0021 

.0001 
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a low  of  7.23  for  the  PVI7  to  a high  of  7.52  for  the  PVI6.  Again, 
the  differences  between  the  yield  relationship  defined  for  the  Vi's 
was  small.  Little  variation  in  the  independent  variables  entering 
the  regression  relationships  was  found  due  to  the  relatively  high 
correlation  between  the  Vi's. 

The  independent  variables  proving  to  be  the  most  significant  at 
harvest  were  the  VI  multiplied  by  surface  soil  moisture  and  sub- 
surface soil  moisture.  Other  variables  proving  to  be  significant 
at  harvest  were  the  code  variable  differentiating  between  wheat 
types  and  the  sum  of  surface  and  subsurface  soil  moisture  raised  to 
the  fourth  power.  The  sum  of  surface  and  subsurface  soil  moisture 
variable  raised  to  the  third  power  was  found  to  be  significant  in 
the  LAI  relationship  while  the  VI  added  to  subsurface  soil  moisture 
was  found  to  be  significant  in  the  TVI6  relationship. 

3.3  RELATIONSHIPS  TO  PLANT  DENSITY  AND  PLANT  HEIQIT 

3.3.1  Plant  Density.  The  relationship  between  plant  density  and  the 

independent  variables  produced  relatively  strong  relationships  at 
tillering  and  sten  extension.  The  analysis  was  limited  to  these 
two  growth  stage  intervals,  because  of  the  method  used  to  record 
the  plant  density  data  (refer  to  Section  1.4.2). 

The  coefficient  of  determination  and  standard  deviation  between 
the  observed  and  estimated  plant  density  at  tillering  were  .62 
and  .46,  respectively  (Table  3-8).  Those  variables  which  were 
significant  in  explaining  the  variation  in  plant  density  at  tiller- 
ing included  the  AVI,  DVI , PVI7,  growth  stage,  surface  and  sub- 
surface soil  moisture.  All  of  the  Vi's  were  included  in  the  list 
of  independent  variables.  Interactive  terms  similar  o those  defined 
for  the  yield  analysis  were  not  defined  for  the  plant  density  or 
plant  height  analysis. 

Growth  stage  was  identified  as  the  most  significant  variable 
followed  by  the  AVI  and  surface  soil  moisture.  The  tillering 
growth  stage  interval  spans  growth  stages  1-5  during  which  sig- 
nificant changes  in  plant  density  occur.  During  tillering  the 
plant  density  is  low  and  the  soils  greatly  affect  the  reflectance 
responses  recorded  by  Landsat.  The  Vi's  cannot  effectively 
measure  the  increase  in  plant  density  due  to  the  resolution  of 
the  Landsat  allowing  the  predominance  of  soil  response  over  plant 
response. 

At  stem  extension  the  coefficient  of  determination  and  the  standard 
deviation  were  .74  and  .48,  respectively  (Table  3-8).  The  variables 
which  were  significant  in  explaining  the  variation  in  plant  density 
included  the  TVI6,  LAI,  wheat  type  and  subsurface  soil  moisture. 

TVI6  was  identified  as  being  the  most  significant  variable  followed 
by  the  code  variable  differentiating  between  wheat  types,  the  LAI 
and  subsurface  soil  moisture.  The  predominance  of  soils  over  plant 
matter  was  negligible  during  stem  extension.  Those  Vi's  computed 
frcm  band  6 seemed  to  measure  plant  density  more  accurately  than 
those  computed  from  band  7.  Both  the  TVI6  and  LAI  are  based,  in 
part,  on  band  6. 
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TABLE  3-8  PLANT  DENSITY  REGRESSION  RELATIONSHIP  AND  ASSOCIATED 
STATISTICS  AT  TJLLERINC  AND  STEM  EXTENSION 


N 

R2 

STD  DEV 

B-VALUE 

STD  ERROR 

F-VALUE 

PROB  >F 

TILLERING 

196 

.62 

.46 

- 

- 

52.31 

.0001 

VARIABLES 

Intercept 

- 

_ 

- 

1.65 

- 

- 

- 

AVI 

- 

- 

- 

0. 14 

0.03 

5.84 

.0001 

DVI 

- 

- 

-0. 19 

0.04 

3.25 

.0001 

PVI7 

- 

- 

- 

0.  20 

0.09 

1.05 

.0278 

CROWSTA 

- 

- 

- 

0.42 

0.  03 

33. 87 

.0001 

SOILMOl 

- 

- 

- 

0.46 

0. 10 

4.60 

.0001 

S01LM02 

- 

-0.08 

0.03 

1.53 

.0080 

STEM  EXTENSION 

31 

.74 

.48 

- 

- 

18.96 

.0001 

VARIABLES 

Intercept 

- 

- 

- 

-39.29 

- 

- 

- 

TVI6 

- 

- 

- 

44.39 

12.54 

12.53 

.0015 

LAI 

- 

- 

- 

-0.20 

0.07 

7.18 

.0126 

Wheat 

- 

- 

- 

-0.58 

0.22 

7.30 

.0120 

S0ILM02 

-0.37 

0.14 

6.79 

.0149 

% 

V 
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3.3.2  Plant  Height.  Relatively  strong  relationships  between  plant 

height  and  the  independent  variables  were  defined  at  tillering 
and  stem  extension.  The  coefficient  of  determination  and  the 
standard  deviation  between  the  observed  and  estimated  plant  height 
at  t i 1 ler ing  were  .66  and  1.40,  respectively  (Table  3-9).  Those 
variables  which  were  significant  in  explaining  the  variation  in 
plant  height  at  tillering  included  the  LAI,  growth  stage,  wheat 
type,  and  surface  soil  rroisutre.  All  of  the  Vi's  were  included  in 
the  list  of  independent  variables.  Interactive  terms  were  not 
defined  in  this  analysis. 

Growth  stage  was  identified  as  the  mast  significant  variable 
followed  by  the  LAI,  the  code  variable  differentiating  between 
wheat  types  and  surface  soil  moisture. 

At  stem  extension  the  importance  of  the  Vi's  in  explaining  the 
variation  in  plant  height  grew  stronger.  The  coefficient  of 
determination  and  standard  deviation  at  stem  extension  were  .91 
and  1.93,  respectively  (Table  3-9).  The  variables  which  proved 
to  be  significant  were  the  PVI7,  TVI6,  growth  stage  and  subsurface 
soil  moisture.  The  PVI7  proved  to  be  the  most  significant 
variable  followed  closely  by  growth  stage. 
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APPENDIX  A 


VEGETATIVE  INDEX  TRANS  FORMATIONS 


AVI  = 

2CH*4  - CH2 

DVI  = 

2.4CH4  - CH2 

GVI  = 

-.283CH1  - . 660CH2  + .577013  + .388CH4 

LAI  = 

4 1 . 325CH1  42  -45.™ 

' CH2  ' CH3 

KVI  = 

GVI  - SOIL  LINE** 

PVI6  = 

(-.498-.457CH2+.498CH3)2+(2>734+>498CH2_>543CH3)2 

PVI7  = 

(.355CH4-.149CH2)2+(.355ai2-.852CH4)2 

TVI6  = 

CH3-CH2  | 1 0 

CH3+CH2 

TVI7  = 

CH4-CH2  0 

CH4+G12 

* CHANNEL  1 = QIl  = BAND  4 
CHANNEL  2 = CH2  = BAND  5 
CHANNEL  3 = CH3  = BAND  6 
CHANNEL  4 = CH4  = BAND  7 

**  The  value  computed  by  the  soil  line  equation  is  equal  to  a constant 
of  "-7"  after  the  Landsat  data  have  been  sun  angle  and  haze  corrected 
The  value  computed  by  the  soil  line  does  vary  without  sun  angle 
and/or  haze  corrections. 
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CORRELATIONS  AND  SIGNIFICANCE  LEVELS  AT  HEADING 
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CORRELATIONS  AND  SIGNIFICANCE  LEVELS  AT  STEM  EXTENSION 
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CORRELATIONS  AND  SIGNIFICANCE  LEVELS  AT  STEM  EXTENSION,  HEADING  AND  FLOWERING 
FOR  FIELDS  GREATER  THAN  OR  EQUAL  TO  30  PIXELS 
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LIST  OF  INDEPENDENT  VARIABLES* 


A 

VI 

(VEGETATIVE  INDEX) 

B 

GRCWSTA 

(GRCWIH  STAGE) 

C 

WHEAT 

(DIMTY  VARIABLE  DIFFERENTIATING  BETWEEN  WHEAT  TYPES) 

D 

SOILMOl 

(SURFACE  SOIL  MOISTURE  IN  INCHES) 

E 

S0ILMD2 

(SUBSURFACE  SOIL  MOISTURE  IN  INCHES) 

F 

SOILM03 

(SUM  OF  SOILMOl  AND  SOILM02) 

G 

SMI  SO 

(SOILMOl  RAISED  TO  THE  SEOOND  PGtfER) 

H 

3T1CU 

(SOILMOl  RAISED  TO  THE  THIRD  POWER) 

I 

SvllFO 

(SOILMOl  RAISED  TO  THE  FOURTH  POWER) 

J 

SM2SQ 

(SOILM02  RAISED  TO  THE  SECOND  PCWER) 

K 

SM2CU 

(SOILM02  RAISED  TO  THE  THIRD  POWER) 

L 

SM2FO 

(SOILM02  RAISED  TO  THE  FOURTH  PQ\ER) 

M 

SM3SQ 

(SOILM03  RAISED  TO  THE  SEOOND  POWER) 

N 

SM3CU 

(SOILM03  RAISED  TO  THE  TITIRD  POWER) 

0 

SM3FO 

(SOILM03  RAISED  TO  THE  FOURTH  POWER) 

P 

VIMM 

(PRODUCT  OF  VI  AND  SOILMOl) 

Q 

VIMSM2 

(PRODUCT  OF  VI  AND  SOILM02) 

R 

VIMSM3 

(PRODUCT  OF  VT  AND  SOILM03) 

S 

VI  AM 

(SUM  OF  VI  AND  SOILMOl) 

T 

VIASM2 

(SUM  OF  VI  AND  SOILM02) 

U 

VI  AM 

(SUM  OF  VI  AND  SOILMD3 ) 

V 

VIM3S0 

(PRODUCT  OF  VI  AND  GRCWSTA) 

* The  independent  variables  listed  were  used  in  the  stepwise  regression 
procedure  corrpleted  for  each  of  the  seven  growth  stage  intervals.  A 
single  vegetative  index  and  its  associated  variables  were  used  as 
inputs  to  the  stepwise  procedure.  The  same  associated  variables  were 
used  in  the  stepwise  procedure  completed  for  each  VI. 


ecp  x 8 1979 

NASA-JSC 


E-l 


PLOT  OF  LAI *GROWST A LEGEND : A = 1 OBS,  R = 2 OBS*  ETC 


v 


I 

I 

I 

oLl  O lo  m nj  rvj  rsj  o < Oj  < + (\j 

» — 

i 

r 

i 

i 

< | 

XXXXO  ►-lDO^-'LuU-  C | 

< <3  <OL_>U_OLjJ*— ibJb^  bJCDCX  <ICJO<  I 

<<<<UUlLliJU^OOUIZrili.^UJ(DOQ<  | 


< 


< < OX  OLl ")  2 UJU.  ' OX O ' 
<<  < < 

< < 

< 

<f  < < <u  CJ 

<<  x oxx<<j«s 


IC2UI 

x cc 


< 


ox  < 


< 


I 

♦ o 


<3 


< 


< < < < < 


< 


<1  <1  < 


< 


< < < <■ 


X < X Lu  X O c < 


I 

I 

I 

f 

« 

t 

I 

I 

4 0s 

I 

I 

i 

I 

I 

I 

I 

I 

I 

♦ cc 

I 

l 

I 

I 

I 

I 

I 

I 

I 

4 

I 

i 

i 

l 

I 

l 

t 

t < 

I k- 

4 sc  m 

I j£ 


i 

\ 

\ 

\ 

I 

4 VJ'I 

i 

\ 

1 

I 


i 


1. 


X <.  < < <L 


< <r  x 


I 

I 

I 

I 

-f 

! 

I 

» 

\ 

I 

1 


< 


FIGURE  2-4  YIELD  VS.- LAI  AT  STEM  EXTENSION 
AI  LEGEND:  A = 1 OHS*  H = ?_  OBS*  ETC 


FIGURE  2-5 YIELD  VS.LAI  AT  HEADING 

PI  OT  OF  YIELD»L'\I  LEGEND:  A = I OBS,  B = ?_  OBS»  ETC 


FIGURE  2-7  YIELD  VS*.  LAI  AT  RIPENING 
PLOT  OF  Y I F L 0 A I LEGEND : A = 1 ORS,  H ? OBS.  ETC 


FIGURE  2-8  YIELD  VS,  LAI  AT  HARVEST 
PLOT  OF  Y ILL  0 L A I LEGEND:  A = I OBS,  R = ?_  ORS.  ETC 


FIGURE  ?.-<J  YIELD  VS.  SUBSURFACE  SOIL  MOISTURE  AT  PLANTING 
PLOT  OF  Y I F.LO^SO  I LMO?  LEGE  NO : A = \ DBS,  B = ? ORS,  ETC. 


FIGURE  2-10  YIELD  VS.  SUBSURFACE  SOIL  MOISTURE  AT  TILLERING 
PLOT  OF  Y IEL0«S0ILM02  LEGEND : A = 1 ORS.  B = P ORS.  ETC. 
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FIGURE  2-12  YIELD  VS  SUBSURFACE  SOIL  MOISTURE  AT  HEADING 
OS0ILM02  LEGEND:  A = i ors,  R = 2 OHS.  ETC. 


FIGURE  2-13  YIELD  VS, SUBSURFACE  SUIL  MOISTURE  AT  FLOWERING 
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NOTE:  15  OBS  HIDDEN 
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